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Abstract: We determine the most general form of the equations of relativistic su- 
perfluid hydrodynamics consistent with Lorentz invariance, the Onsager principle and 
the second law of thermodynamics at first order in the derivative expansion. Once 
parity is violated, either because the U(l) symmetry is anomalous or as a consequence 
of a different parity-breaking mechanism, our results deviate from the standard text- 
book analysis of superfluids. Our general equations require the specification of twenty 
parameters (such as the viscosity and conductivity). In the limit of small relative su- 
perfluid velocities we find a seven parameter set of equations. In the same limit, we 
have used the AdS/CFT correspondence to compute the parity odd contributions to 
the superfluid equations of motion for a generic holographic model and have verified 
that our results are consistent. 
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1. Introduction 

In this paper we work out a complete theory of relativistic superfluid hydrodynamics 
including dissipative terms at first order in a gradient expansion, allowing for parity 
violation and the presence of triangle anomalies. In other words we work out the most 
general form of the equations of 3+1 dimensional s-wave superfluid hydrodynamics con- 
sistent with Lorentz invariance and the second law of thermodynamics. While we work 
in a relativistic context throughout this paper, our final results admit a straightforward 
non relativistic limit, and are easily formulated in a non relativistic context. 

The theory of superfluid hydrodynamics has a long history. The equations of ideal 
superfluidity (i.e. superfluid dynamics in the absence of dissipative terms) were worked 
out over 60 years ago by Landau and Tisza [|I|, Q in a non relativistic setting. They were 
generalized to a relativistic superfluid in the early 80 's by Israel and Khalatnikov and 
Lebedev || |], || and later reformulated by Carter and Khalatnikov || [7| and by Son || . 
In most of this work we will use a formulation of superfluid dynamics close to that of || . 
In a beautiful recent work, Sonner and Withers (see also [|1(|) have used the equations 



of Einstein gravity to rederive the Landau-Tisza equations for superfluids that admit a 
dual gravitational description via the AdS/CFT correspondence of string theory. This 
development yields independent evidence for the correctness and completeness of the 
Landau-Tisza theory of ideal superfluidity. 

The focus of the current paper is on the one derivative (dissipative) corrections 
to the Landau-Tisza equations. Dissipative corrections to relativistic superfluids have 
been extensively studied in the literature [|ll], 12[ [13, [L4|, 15, [16|, [TUJ and they generalize 



(and extend) the textbook derivation of such corrections in the non-relativistic limit |17 



|19|j. While it is straightforward to list the most general dissipative corrections to the 
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stress tensor, charged current and Josephson relations allowed by Lorentz invariance, 
it turns out that such a listing also allows for many unphysical possibilities. It is a 
physical requirement that any hydrodynamic flow be equipped with an entropy current 
Jii. The second law of thermodynamics requires that the increase in the entropy in any 
compact, spacelike, region be larger than the incoming entropy flux through the surface 
of that region; this requires that at every point in spacetime and for every conceivable 
fluid flow 

<V£>o. (l.i) 

An interesting and important fact about fluid dynamics, whether superfluid or not, is 
that the requirement that the divergence of the entropy current always be positive gives 
rise to important constraints on dissipative corrections to the equations of motion. 

As far as we are aware, all previous studies of dissipative corrections to the equa- 
tions of superfluid hydrodynamics assume, on intuitive grounds, that the entropy cur- 
rent of superfluid hydrodynamics takes a particular canonical form, Jq . These 

±r J J i ' o canon 

studies then determine the dissipative corrections to the energy momentum tensor and 
charged current (and Josephson condition) consistent with the positivity of the diver- 
gence of this canonical entropy current and with covariance under time reversal. The 
latter restriction is usually called the Onsager relations. 

In [201, Son and Surowka observed that, in the presence of triangle anomalies, the 



entropy current of ordinary fluids (non superfluids) deviates from its canonical form. 
This observation makes clear that the intuition that fixes the entropy current to its 
canonical form is not infallible. For this reason, the starting point of our analysis in 
this paper is the assumption that an entropy current with positive definite divergence 
exists. However, we make no assumption about the form of this current beyond the 
requirements of symmetry. 

More formally, we allow the (postulated positive divergence) entropy current to 
take the form 

= + (1-2) 

i 

where Jo is the canonical entropy current referred to above, is a basis of on- 

J canon ± */ > % 

shell inequivalent one derivative vectors and Vi are (initially) unconstrained coefficient 
functions. We then demand that the divergence of J$ be positive semi-definite for 
any solution of the equations of superfluid hydrodynamics on an arbitrary background 
spacetime and that the Onsager relations are satisfied. 1 These restrictions fix most of 



1 Recall that the second law of thermodynamics must apply in any conceivable consistent situation. 
In particular it must apply when the system is formulated on an arbitrary background spacetime 
provided the system is free of diffeomorphism anomalies. This condition is true of all experimental 
superfluids as well as all superfluids obtained via the AdS/CFT correspondence. 
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the ViS in ( |1.2|) and also restrict the possible transport coefficients of the theory. 

As we discuss in detail throughout this work, the requirement that the entropy 
current be of positive divergence in an arbitrary background spacetime provides powerful 
constraints on the form of the entropy current and through it, on the form of the possible 
dissipative corrections to the hydrodynamic constitutive relations even in flat space. For 
instance, the divergence of the entropy current could contain a term proportional to 

V M J£ oc ViR^uW + ... (1.3) 

where is the fluid velocity, R^ v is the Ricci tensor, and V\ is one of the coefficient 
functions in ( |1 . 2|) . The divergence of the entropy current may also contain many other 
terms proportional to v,- t and independent of curvatures. However, for any given fluid 
flow these other terms can be held fixed while R^u^u" is made arbitrarily negative by 
tuning the curvature tensor. 2 It follows that the divergence of the entropy current is 
positive for an arbitrary fluid flow on an arbitrary spacetime only if v\ = 0. Thus, we 
find a constraint on the entropy current for fluid motion in a flat space background, 
even though we needed to move to a curved spacetime in order to obtain this constraint. 

The result of our extended analysis is as follows. For superfluids which are parity 
invariant we find that up to a physically unimportant ambiguity (see sections ||] and |] 
for details) all the v^s in (|1.2j) should be set to zero — the entropy current agrees with the 
canonical entropy current. 3 This result provides a check of the intuition reviewed in, say, 



17 , that the entropy current should take its canonical form. It follows that the most 



general structure of dissipative terms in parity preserving superfluids is given by the 



14 parameter family described in [10] which generalized the 13 parameter construction 
of Clark and Putterman []I8 , [19j ) . 

Let us now turn to the case of superfluids that do not respect parity symmetry. In 
this case we find that the entropy current is not constrained to take the canonical form. 
Not only are the ViS non zero, two of them remain undetermined (meaning that they 
are completely free functions of local thermodynamical variables). In additional there 
are two physically unimportant ambiguities which we describe in detail in sections |3| 
and g. 

It turns out that an 18 parameter set of equations of superfluid hydrodynamics 
are consistent with the positivity of any given choice of entropy current. 14 of these 
parameters multiply parity invariant structures; these are the parameters that pa- 
rameterize superfluid hydrodynamics in parity preserving systems. The remaining 4 



2 Note that curvature tensors do not contribute to the fluid equations at first order, so it is consistent 
to hold fluid flows fixed while taking curvatures to be very large. 

3 As a warm up for this analysis we show in section ^] that a similar result is true of parity invariant 
charged fluid dynamics in the absence of superfluidity 
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parameters multiply parity odd structures, and are completely new. Thus the most 
general equations of parity odd superfluid hydrodynamics are parameterized by 14 +4 
+2 parameters; the last two parameters being the undetermined f^s in the entropy 
current. 

At this point it would be useful to clarify some of the terminology we use. We have 
stated that the equations of motion of superfluid dynamics require the specification of 
twenty unknown functions of the thermodynamic variables. This may be contrasted 
with the equations of motion of a normal, charged conformal fluid which requires the 
specification of two parameters, the shear viscosity and conductivity. For the normal 
charged conformal fluid, the shear viscosity and conductivity also control the amount 
of entropy produced by the fluid and so, these parameters are also called dissipative 
parameters. In the case of parity violating superfluids, five of the twenty parameters 
are not associated with entropy production so, only fifteen of the twenty parameters 
are dissipative in the sense described above. In normal charged conformal fluids the 
conductivity controls not only the response of the current to an external electric field, 
but also the response of the current to changes in temperature and changes in chemical 
potential. The situation in superfluids is similar: the twenty parameters, and in par- 
ticular the two undetermined parameters associated with the entropy current, control 
the response of the system to various changes in the hydrodynamic variables. 

When the superfluid velocity is too large superfluidity breaks down. For this reason, 
when considering experimentally accessible superfluids such as liquid Helium it is often 
particularly interesting to study dissipative corrections in the limit that the normal and 
superfluid velocities are collinear. In this limit the 14 parameter family of dissipative 
corrections of parity invariant superfluids reduce to the five parameter family described, 
for instance, in the classic text book of Landau and Lifshitz |17|. 4 When parity is 
violated, the twenty parameter family of solutions reduces to seven. The two additional 
parameters determine (in an appropriate frame) the response of the charged current to 
a magnetic field and to "chiral vorticity", i.e., changes in u^ 1 = \e^ va ^u v d a up. 

When the superfluid density is set to zero we are in the normal phase. In this limit 
the two undetermined parameters of the superfluid theory reduce to an integration 
constant and the contributions to the charge current proportional to the magnetic field 
and chiral vorticity are completely determined by the triangle anomaly POI (up to an 
integration constant). In such a limit, we recover the results of [^T], ^2|, ^(J (see also 

We emphasize that, as opposed to the configuration at zero superfluid density, the 6 
new arbitrary parameters that appear in parity non invariant superfluid hydrodynamics 



4 These reduce further to only three parameters in the case of a conformally invariant theory. 
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may be non vanishing even in the absence of an anomaly. This difference may be of 
significance; practically speaking, anomalies are intrinsically relativistic phenomena 
and always vanish in a nonrelativistic setting. Consequently, the results of this paper 
suggest that a term in the particle number current proportional to the vorticity could 
possibly show up in non relativistic table top experiments involving superfluids or 
superconductors which violate parity (e.g. non centro symmetric superconductors). 

We were led to the study of dissipative effects in superfluids in order to understand 
the results of certain holographic computations using the AdS/CFT correspondence of 
string theory. We emphasize, however, that our eventual derivation of the general equa- 
tions of superfluid hydrodynamics makes no use of the AdS/CFT correspondence and 
so, applies to all superfluids, not just those that admit a dual gravitational description. 

In the second part of our paper we test the collinear limit of our general formalism 
using the AdS/CFT correspondence. In the case of parity preserving superfluids, the 
authors of |16|, [UJ showed that it is possible to use the fluid gravity map [FJ|, ^5], |26], |27|, 
[28| , |2T| , p2| , to derive dissipative corrections to the equations of superfluid dynamics 
by extending the analytically tractable superfluid model of PU| which is valid close to 
the phase transition point. The constitutive relations so obtained fit within the 14 
parameter fluid dynamical framework spelt out in [T(J which generalize the text book 
predictions of |18], 0)- 111 particular, it was explicitly verified that the entropy current 
obtained by holographic methods matches the canonical entropy current in agreement 
with the intuition of []i~7| . 

In this paper we consider a generic asymptotically AdS$ gravitational system which 
involves a Chern Simons term for a single bulk U(l) field. Such a Chern Simons term 
signifies the presence of a U(l) 3 triangle anomaly in the dual field theory and so, in 
particular, introduces parity violating terms into the effective superfluid dynamical 
description. Particular truncations of type IIB supergravity for which our analysis is 



valid can be found in [[31], |32fl. We find that, in the collinear limit, it is possible to 
obtain integral expressions for all the parity odd transport coefficients in terms of the 
background solution. Our results are in perfect agreement with the collinear limit of 
the modified theory of parity violating superfluid dynamics which we present in section 
[|. We regard this agreement as a nontrivial (though as yet limited) check of the general 
results of section f|. It should be possible to use the AdS/CFT correspondence away 
from the collinear limit to find much more extensive test our results; we leave this to 
future work. 

Note Added: After the work presented in this paper was completed, we received 
a paper |33|] that has substantial overlap with section | of this paper. While the 
approach of |3l| is similar to the one adopted in this paper, our final results differ in 



several qualitative and quantitative respects. We present a brief comparison with [33 
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in section ||. 



2. The theory of charged fluid dynamics 

In this pedagogical section we construct the most general equations of Lorentz invariant 
charged fluid dynamics consistent with the second law of thermodynamics. Our goal 
is to illustrate our method for determining the most general form of fluid-dynamical 
equations of motion in a simple and familiar context before tackling the slightly more 
complicated case of superfluids. The final results of this section are well known; the 
novelty of this section lies in our method of computation. 

The long-wavelength degrees of freedom of a locally equilibrated system with a 
single global U(l) charge can be taken to be the velocity field u^(x) (normalized so 
that u^u^ = —1), the temperature field T{x) and a chemical potential field Both 
the energy momentum tensor and the charged current can be written in terms of these 
five fields and their gradients. The equations of motion of charged fluid dynamics are 
the conservation of the stress tensor and charge current 



which provides the five equations for the five hydrodynamic fields. In these equations 
we have allowed for the possibility that the current in question has a U(l) 3 anomaly. 
We call the coefficient c the anomaly coefficient. We have also allowed the current to 
be coupled to an external source with field strength F^ v . To completely determine the 
equations of motion it remains to determine the dependence of T^ v and on the fields 
u^(x), T(x), fi(x) and their derivatives. 

By considering a stationary fluid for which = (1, 0, 0, 0) and using boost in- 
variance one can argue that the stress tensor and charge current take the form 



where T^ ss and J^ iss are the contributions to the stress tensor and charge current that 
involve derivatives of /x, T and u^. The equations that express T£? ss and J% iss in terms 
of fluid dynamical fields and their derivatives are termed constitutive relations. In the 
long wavelength fluid dynamical limit it is sensible to expand the constitutive relations 
in powers of derivatives of the fluid dynamical fields w M , T and /j. We will refer to 
such an expansion as a derivative expansion and refer to the terms which are linear 
in gradients as first order terms. In this paper we work only to first order in the 



(2.1) 



= (p + P)u»u" + PtT + T£ 



(2.2) 
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derivative expansion. The electromagnetic source term F^ v is taken to be of first order 
in derivatives in this counting. 



Field Redefinitions and frame choices 

Note that the fluid temperature T, chemical potential \i and velocity are ther- 
modynamical concepts that are well defined in equilibrium but have no microscopic 
definitions in dynamical situations. In other words, we are always free to redefine the 
thermodynamic variables into primed ones according to the equations 

= u'v + 5u» 

T = T' + 5T (2.3) 
H — fx' + 5/J, 

where is an arbitrary one derivative vector that obeys Su^u^ = 5u^u'^ = and 
5T and 5fi are arbitrary one derivative scalars. The primed and unprimed fields are 
each equally good definitions of the velocity, temperature and chemical potential fields. 
Physically meaningful assertions, such as the constitutive relations for T^ ss and J^ ias , 
must only involve field redefinition invariant quantities. 

Thus, let us determine the field redefinition invariant combinations of T^ ss and 



J diss- Under the field redefinition ( [2. 3D 



STZs = {^5u v + u W) (P + P ) + u Vd(P + p)+ rTdP 
ZJLs = ^ + dqu» 



where 



iH v 
diss 


rpl PlV 

diss 


rpjiV 
diss 


diss 


_ j/fj, 
diss 


- 

diss 



and df(fi,T) represents the change in the function / under the first order variable 
change ( |2.3| ). It useful to decompose Tjj? s and J% iss into 5*0(3) invariant tensors, vectors 
and scalars The S0(3) that we are referring to is the group of rotations orthogonal to 
m m . To this end we introduce the projection matrix 

P i*> = v »» + u » u v . (2.5) 
We find that there is one tensor, two vectors and three scalars. The unique tensor 



p" a p u A - -g-JWfSL (2-6) 
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is automatically field redefinition invariant. The two vectors P^ a Pdiss u P an d P^ a J a 



(2.7) 



transform under field redefinitions as 

S (P^Lup) = -{P + p)8u» 
5 (P^r) = q5u» 

so that the unique invariant combination of vectors is given by 

P"J a + (P^lup) • (2.8) 

The three scalars transform under field redefinitions as 

<* I PafitfZ.) = 3dP 



6 ( u <*Tdf ss u^ = dp ( 2 - 9 ) 
S (u a J a ) = —dq 

so that the unique invariant scalar is given by 

\ (P*pT2£.) - ^ (u a Tf ss u,) + ^ {u a r) (2.10) 

where ^ is taken at constant q and |p is taken at constant p. 

Instead of working in a manifest field redefinition invariant manner, it is sometimes 
convenient to 'fix' the field redefinition ambiguity by imposing five additional conditions 
on the thermodynamic fields so that they are well defined. Different choices of fixing the 
ambiguity are referred to as frames. One often used frame is the so called Landau frame, 
in which the velocity and temperature fields are defined to obey the conditions T^ ss u v = 
and J^ iss u u = 0. This gives one vector and two scalar conditions, matching the field 
redefinition degrees of freedom. Another choice of frame is the Eckart frame which is 
defined by the conditions J% iss = together with u^T^ ss u u = 0. The expressions for 
the invariant vector ( f2.8| ) and the invariant scalar ( 2.10|) greatly simplify in either of 



these frames. In this paper we adopt no such 'gauge' choice but work in a fully field 
redefinition invariant manner. 



The strategy for the rest of this section 

In order to complete our specification of the equations of charged fluid dynamics, we 
need to specify T^ ss and J% iss (or more precisely the field redefinition invariant parts 
Q2.6Q , ( p.8|) and (|2.10|) of these expressions) as a function of first derivatives of fluid 
dynamical fields. Of course, in any particular dynamical system, the explicit form 
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of the constitutive relations for T^ ss and J^ iss can be determined only by a detailed 
dynamical computation. In this paper we will be interested not in computing the 
precise form of these quantities in any particular system, but in parameterizing the 
most general form that the constitutive relationship can take in any system. As we 
will see below, it will prove possible to completely determine the form of the first order 
constitutive relations up to three undetermined dissipative parameters, each of which 
is an arbitrary function of T and p. 

We proceed as follows. As in any effective field theory, we start by writing down 
all possible expressions which may contribute to T^ ss and J^ isa - We then eliminate 
those that do not satisfy the symmetries of the theory, Lorentz invariance in this case. 
In addition, since we are dealing with a hydrodynamic theory, we must ensure that 
the second law of thermodynamics is satisfied. As explained in the introduction, we 
demand the existence of an entropy current of positive semi-definite divergence even 
when the theory is formulated on a curved background. As alluded to in section [1], the 
entropy current is defined to be a four vector J$ satisfying two requirements. The first 
is that in a configuration where the fluid is in uniform motion, 5 

Jg = su 11 (for a spacetime independent configuration) (2-11) 

with s the entropy density which is related to p, P, q, p and T through 

p + P = sT + pq. (2.12) 

Our second requirement of the entropy current is that its divergence is positive semi- 
definite in an arbitrary curved background, 

V^J£ > (2.13) 

implying that the entropy increase in any region is always greater than the entropy 
inflow into that region. 

For a perfect fluid level (i.e. a fluid in which all gradient terms have been neglected — 
Tdiss = ^ diss = 0) the entropy current is given by fl2.11| ). At this order it is not difficult 



to verify that V M (su^) = using (pT2f) and dP = sdT + qdp. 

Once the gradients of v?, T and p/T are non vanishing the divergence of the entropy 
current no longer vanishes. Indeed, the divergence of the entropy current at the one 
derivative level will be the focus of much of the rest of this paper. We will demand, on 



5 Such a configuration is a stationary, dissipation free solution to the equations of fluid dynamics. 
Indeed it may be obtained by boosting a uniform fluid at rest (by which we mean a uniform fluid with 
velocity field u» = (1, 0, 0, 0)). 



- 10 - 



physical grounds, that it is possible to modify (|2.11|) by first order corrections so that 
( |2.13|) will be satisfied. This requirement will turn out to constrain the possible forms 



of T dL and J Ls- 



We start our analysis in section |2J] by considering parity conserving charged fluids. 
In section |2.2| we move on to describe parity- violating fluid dynamics. 



2.1 Parity Invariant Charged Fluid Dynamics 

Consider a hydrodynamic theory in the presence of external electromagnetic fields sat- 
isfying with c = 0. Following the general prescription described at the beginning 
of this section, we would like to write the most general parity-invariant and Lorentz 
invariant contributions to J^ iss , T^ ss and Jg which involves a single derivative of the 
hydrodynamic fields u^, T and [i. This is carried out in section |2.1.1| . We then work out 
the the restrictions on these terms by requiring that the entropy current has positive 



semi-definite divergence. This is described in section 2.1.2 



2.1.1 Classification of one and two derivative data 

We begin our analysis on a technical point. The tangent space about any point in our 
spacetime manifold has an SO(3, 1) rotational invariance. However, the fluid velocity 
vector, u^(x), takes a definite value at that point and breaks this rotational group down 
to SO (3). It is useful to decompose all derivatives of fluid dynamical fields, at any given 
spacetime point, into representations of this residual SO(3) rotational group. 

In the first column of table [l] we have classified all expressions formed from a single 
derivative of any of T(x) and fi(x) according to their SO(3) and parity transfor- 

mation properties. We refer to these expressions as one derivative fluid dynamical data. 
We have also classified one derivative expressions constructed out of the background 
electromagnetic fields according to their SO (3) and parity transformation properties. 
We will refer to these as background data. As fluid and background field data enter 
our analysis on an even footing, we have listed these expressions together in the first 
column of table [TJ. 6 

Not all the expressions in the first column of table |I] are independent under the 
equations of motion. The equations of motion can be used to solve for some pieces of 
data in terms of other data. The classification of the equations of motion according to 
their 5*0(3) and parity transformation properties can be found in the middle column 
of table |l|. Note that there are no tensor equations of motion. 

6 Since all curvature invariants built out of the background metric have at least two derivatives, 
there is no one derivative data associated with the metric. 
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In the last column of table p] we have listed a choice of independent data. By this 
we mean a choice of independent one derivative fluid dynamical expressions and one 
derivative field expressions in terms of which all others can be solved for. 



SO (3) and P 
classification 


All data 


Equations of motion 


Independent data 


Scalars 




u^ v T^ = 
V M J" = 


Si = 


Vectors 


P» v d v T 
F^u v 


P^d p T pv = 


v x = p^d^+ F ^ 

V 3 = F^uu 


Tensors 








Pseudo vectors 


\& vafi u v d a up 




w m = l e ^u v d a u p 

B v = l_ e ^ uyFaf} 



Table 1: One-derivative expressions classified according to their transformation laws under 
the SO(3) residual symmetry and parity. The first column lists all one derivative data. The 
second column lists the equations of motion. The last column lists a choice of independent 
data. See (2.14) and ( |2.15| ) for the definition of and P pu respectively. 



While some of the expressions used in table |l| such as 

<V = X -P^P^ (V a up + V p u a - P a p {V x u x )) 



and 



(2.14) 



(2.15) 



are standard, some of our notation isn't. The new notation has been introduced in 
order to prepare the reader for later sections. In particular V3 is the electric field in the 
rest frame of a fluid element. In the conventions of Son and Surowka lEDI we have 



V, = E„ 



(2.16) 



We will soon construct an entropy current that includes terms which are first 
order in derivatives. The divergence of such an entropy current is of second order in 
derivatives and includes terms quadratic in first order fluid (and background field) data 
plus expressions built out of two derivatives acting on fluid fields or single derivatives 



- 12 - 



of electromagnetic field strengths. We refer to the second class of expressions as two 
derivative scalar data. When studying the divergence of the entropy current it is useful 
to have a listing of independent scalar two derivative data. 

In the first column of table |2] we list the most general fluid and background field 
(but not curvature related) two derivative data that transforms as an SO (3) scalar. 
More explicitly, we list all scalar expressions formed by acting with two derivatives 
on w M (x), T(x) and together with all scalars formed from the action of a single 
derivative on electromagnetic field strengths. 7 In the second column of the same table 
we list all scalar two derivative equations of motion. In the last column of the same 
table we list our choice of independent two derivative scalar data (in terms of which we 
have solved for all the other two derivative scalars). 



All data 


Equations of motion 


Independent data 


u»u u d^d u T 
v?u v d l & v (£) 

P^dpdy (f ) 

d^F^u v ) 


V^ V T^ = 
P^V P V^T^ = 
u u V u V^J^ = 





Table 2: Parity even two derivative scalar data for charged fluids. The first column lists all 
six second order scalars constructed from two derivatives of the hydrodynamic variables and 
background field strengths. The second column lists the three scalar two derivative equations 
of motion. The last column lists one choice of a 6— 3 = 3 dimensional basis for the independent 
two derivative scalar data. 



2.1.2 The general entropy current and its divergence 

Armed with the listings in tables [I] and |2] we now proceed with our analysis. Traditional 
studies of first order charged fluid dynamics (see, for example, flT7| ) assume that the 
entropy current takes a canonical form, 8 

canon = SV ^ ~ 7^ U ^dlSS ~ T^^dlSS " (2-17) 



7 It is also easy to list two derivative fluid data in the 3, 5 and 7 dimensional representations of 
SO(3), but that will not be required in what follows, so we do not present such a listing. 



8 As explained in JlQ ] the expression in (2.17) is frame invariant, i.e. invariant under a first order 
field redefinition of T, and fi. Note that the second term on the right hand side vanishes in the 
Landau frame while the third term vanishes in the Eckart frame. 
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As we explained in the introduction, in this work we will not make any prior assumption 
about the form of the entropy current. According to the analysis of section [2.1. 1| the 
most general parity even first order entropy current is given by 

3 

^ = ^ anon + *iW + 5>^ (2.18) 

i=i 

where Si and Vi are defined in the last column of table [l], and s\ and the vis are 
arbitrary functions of ^ and T. 

We now explore the constraints obtained by enforcing the positivity of the diver- 
gence of the entropy current ( |2.18|) . It is easily demonstrated (see, for instance, |L7|, fLOl ) 



that the divergence of the canonical part of the entropy current is given by 

V Vscanon = ^/l ^ jT^ ^dlSS ~ ^J^) ^Tp ^ ^diss' (2-19) 



The right hand side of ( [2.19D is a quadratic form in one derivative fluid and background 



electromagnetic field data. The divergence of the non canonical part of the entropy 
current in ( |2.18|) is also a two derivative expression but is composed of two kinds of 



terms. The first set of terms are linear in independent two derivative and curvature 
data. Such terms are always inconsistent with the positivity of the entropy current, 
and so we must choose Si and i>j so that these terms vanish. The second set of terms 
contains products of one derivative terms. Such terms would modify the quadratic form 
on the right hand side of ( p,19[ ) and do not necessarily vanish. Schematically, we have 

q jii. _ ( independent two \ . ( quadratic form m\ ^ 20) 

^ s yderivative and curvature data/ \ first order data / \ ■ J 



The first term on the right hand side of ( |2.20| ) must vanish while the second term must 
be tuned to be positive. 

2.1.3 Constraints from positivity of the divergence of the entropy current 

We will first explore the constraints that follows from the requirement that no two 
derivative data appears in the divergence of the entropy current. As explained pre- 
viously, this implies that the first set of terms on the right hand side of (|2.20|) must 
vanish. We will implement this condition separately for two derivative and curvature 
terms. 

Constraints from the vanishing of 2 derivative terms 

The two derivative part of the divergence of the entropy is given by 

-vxP^V^ + (si + ^KV M cV + (v 3 + ^) V^F^Uu) . 
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This expression is a linear combination of the three independent two derivative pieces 
of data (see Table |2|). It follows that the vanishing of two derivative terms requires 
us to set the coefficients of each of these terms to zero, i.e. to set V\ — v 3 — and 
t>2 = — S\. Thus the vanishing of two derivative terms in the divergence of the entropy 



current restricts the entropy current (2.18) to take the form 



J§ = + «i ( W - Vf) • (2.21) 

where s\ is still an arbitrary function of T and \l. 

Constraints from vanishing of curvature terms 

According to Q2.21 ) the entropy current has a one parameter ambiguity, s\. Were we to 



restrict our attention to a flat space background we would not have been able resolve this 
ambiguity. Consider a charged fluid propagating on an arbitrary curved background. 
The cancellation of two derivative terms proportional to S\ is now incomplete; it is not 
difficult to check that there is an additional, curvature dependent term in the divergence 
of the entropy current proportional to s\R a fiU a vP with R a p the Ricci tensor. This term 
is inconsistent with positivity of the divergence of Jg. Thus, we are forced to set si = 0. 

We conclude that the requirement that the divergence of the entropy current is pos- 
itive in an arbitrary curved background forces the entropy current to take its canonical 



form, justifying the assumptions of standard treatments of fluid dynamics e.g. [|17 



2.1.4 Constraints on dissipative terms 

We have demonstrated that the entropy current takes its canonical form and conse- 



quently that its divergence is given by ([2.191) . It is now not difficult to work out the 



constraints on dissipative terms that ensure the positivity of the quadratic form on the 
right hand side of ( 2.19|) . We outline the calculation here. 



Consider the expansion of V M (%■) and —V\ = — £ which appear on the right 
hand side of ( |2.19| ) into SO (3) invariant tensors vectors and scalars. We find a single 
tensor, a^, two vectors, 

Vf, and (P^^jf- + (u.d)vP 

(see table [I] for a definition of the vector V±) and three scalars, 

(u-d)T 



T 



(u-d)u, and (V-m) 



While the two vectors are completely distinct off-shell, it turns out that the equations of 
motion imply that they are proportional to each other on-shell. Similarly, the equations 
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of motion imply that the three scalars are also proportional to each other on-shell. As 
we demonstrate in Appendix [B] the explicit relations are 



(u.d)T 
T 

{u.d)v = — 



T 



dP 

dp 
dP 
dq 



(V.Ti) 



J P 



d T 

P^^r- + {u.8)u» 



(V.«) 
qT 



(2.22) 



t v ' p+P 

Plugging these relations into ( [2.19 ), we can rewrite the divergence of the entropy current 
in the form 



V M J S 



T 



diss lab 



pab 



d diss * 



3 

q 



dP v dP „ 

-Q^{Ufj,u u T diss ) + -Q^-(u^J diss ) 



( U v T diss) 



l diss°^v 



(2.23) 



P + P 

where Vf*, and were defined in Table [I] and (|2.16|) . We collect their definitions 
here for convenience: 



E„ 



B„ 



Vi, 



1 



' pa/3 



E 



T 



We will now use ( |2.23|) to constrain the constitutive relations of fluid dynamics, i.e. 



the expressions for T^ ss and J diss as a linear expansion in first order scalars, vectors 
and tensors. To first order in gradients there is only one independent scalar data so the 
scalar parts of T^ ss and J diss are necessarily proportional to V -u. The vector parts of 
T%? Ba and J diss must each be expanded as a linear sum of the three independent vectors 
listed in Table [TJ The tensor in Table [I| is proportional to a^ u since there is only one 
50(3) invariant tensor. It follows from group theory that positivity of the divergence 
of the entropy current implies positivity of the scalar, vector and tensor components 
separately. Thus, we have 



piiv 

p/t -pvrpaf3 p 

r a r B 1 diss <- r < 



rpa/3 



a/3 1 diss 



ft i' 



a 



T a 4- 

° diss ' 



■ diss 



diss lab 



P 



ab 



P + P 

dP , „ UU s BP . u 

U ^ T diss) + -HZ^Jdi! 



dp 



Oq 



—rja 

-f3d a u a 



(2.24) 



- 16 - 



where 

77 > 0, k > 0, p > 0. 

These three coefficients are the shear viscosity, 77, the heat conductivity, k, and the 
bulk viscosity, /3. The bulk viscosity is traditionally denoted by ( but in this work we 
reserve ( for different use. 9 



Several aspects of ( |2.24j ) deserve comment. First, the requirement of positivity 
does not individually constrain the three scalar and two vector pieces in and JjL,, 



but only constrains the combinations that appear in ( |2.10| ) and ( |2.8| ). This is exactly as 



we would expect: only field redefinition invariant data can be constrained in a physical 
way. The vectors and scalars that are left undetermined are unphysical; they can 
be changed, or chosen arbitrarily, by a field redefinition. Despite appearances, (|2.24|) 
constitutes a complete determination of the constitutive relations of our system. 

We also note that we could have used the fact that the divergence of the entropy 
current is frame invariant (see to determine the frame invariant scalar, vector and 
tensor combinations in ( |2.6| ), ( |2.8| ) and ( |2.10| ); the expression on the right hand side of 



( |2.23| ) must arrange itself into such frame invariant combinations. 



The third aspect to note is that the constraints of positivity are relatively mild in 
the scalar and tensor sector. The expansion of scalars and tensors is the most general 
one permitted by symmetry; the requirement of positivity merely imposes inequalities in 
the coefficients of this expansion. However, the constraint on vectors is much stronger. 
Symmetry alone would have allowed the expansion of the second line in ( |2.24| ) as an 



arbitrary linear combination of the 3 vectors Vi, V2 and V3. However the requirement 
of positivity sets the coefficients V2 and V3 to zero, 10 apart from imposing an inequality 
on the coefficient of the third. We will see this pattern repeated and magnified in the 
study of superfluid dynamics in sections || and |] in the scalar, vector and tensor sector. 

2.2 Parity non invariant charged fluid dynamics 

Let us now turn to the dynamics of fluids that are not invariant under parity transfor- 
mations. According to table [I] we should allow the entropy current to depend on an 
additional arbitrary pseudo vector. Thus, the most general entropy current for such a 







Note that the speed of sound, c s , is related to the variation of the pressure with respect to energy 
density through ^£ = c 2 s . Using dimensional analysis one can conclude that ^ = in a scale invariant 



theory. It then becomes clear that in a conformal theory the left hand side of the last equality in ( 2.24 ) 
vanishes as it should. 

10 Thc origin of this constraint is the observation that the quadratic form ax 2 + bxy + cxz is positive 
only when b = c = and a > 0. The role of x is played by the vector V , while the roles of y and z are 
played by the other two vectors 
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fluid takes the form 



i=l 

In the parity even sector the divergence of this entropy current is identical to the 
one discussed in subsection IO; the arguments in pTTI go through unchanged and in 



particular the cancellation of two derivative and scalar terms set S\ — Vi — 0. In the 
parity odd sector the divergence of the entropy current receives contributions involving 
the dot product of the pseudo vectors cu M and B^ with ordinary vectors. Positivity 
of the divergence of the entropy current implies that such products vanish. 11 This 



restriction was analyzed in detail by Son and Surowka [^0| who found that it leads to 



p/i purpa/3 ■* p rpa/3 _ „_|«> 



p a (jLs + J^p(n v TZ s )j = «V" + + R B B» (2.26) 

(Tdiss) abP ah 9P 8P 



where 



p 3 

<7u = c— + Tpk 2 + T 2 ki 

p 2 T 
aB = C 2T + 2 h 



2 2 q 3 \ 2 2q \ 2q 
fi - 77 -p + T 1 -fi k 2 -h 



(2.27) 

3p + P n J \~ p + P^J"" p + P'" 1 

( 1 g A T 2 q 

K B = c\fi-- —fi )-— —p-ki 

\ 2p n + P J 2 p + P 

and k\ and k 2 are integration constants. We will now argue that the requirement of 
CPT invariance forces k 2 to vanish. 12 The argument goes as follows. Consider the 
CPT transformation x M — > — x M q — > —q (and so fi — > —ft). Under this transformation 
TZs T dL and -J% iss . Also u " so that w " -w" and -». 

B^. Thus under a CPT transformation it must be that while k b — »■ — /<b- 

Consistency of this requirement with ( |2.27 ) sets k\ = 0. Nothing in our argument 



n We will see later that products of vectors and pseudo vectors do not necessarily need to vanish in 
the case of superfluid dynamics. In the current setup vanishing of such bilinear terms follows from the 
fact that the divergence has no squares of pseudo vectors and contains only a single squared vector. 

12 We thank D. Son for pointing this out to us. 
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requires that k 2 vanish (although it would be intersting to find a specific system with 
k 2 7^ 0; k 2 vanishes in all AdS/CFT computations performed so far). 



of an anomaly forces the entropy current to depart from the canonical form (i.e. Ob and 
are never zero if c is nonzero). Second, it induces new terms in the vector part of the 
constitutive relations, proportional to the vorticity and the magnetic field. Third, the 
new contributions to both the entropy current and the vector part of the constitutive 
relations are completely determined (up to an integration constant that is independent 
of T and \i) by the anomaly. In other words, although the constitutive relations take a 
different form from the parity even case, this change in form is completely determined by 
the anomaly, and we have no new free parameters apart from the integration constant 



3. Parity invariant Superfluid hydrodynamics 

By definition, a superfluid is a fluid phase of a system with a spontaneously broken 
global symmetry. When discussing superfluids this forces us to consider the gradient of 
the Goldstone boson as an extra hydrodynamical degrees of freedom in addition to the 
standard variables w M , T and fi. More precisely, if we denote the Goldstone Boson by 
ip (ip is the phase of the condensate of the charged scalar operator) and we also wish 
turn on a background gauge field then 



represents the covariant derivative of the Goldstone Boson and is an extra hydrody- 
namic degree of freedom. 13 According to the Landau-Tisza two fluid model the super- 
fluid should be thought of as a two component fluid: a condensed component and a 
non condensed or normal component. The velocity field of the normal fluid is given by 
and the velocity of the condensed phase is proportional to £ M . It is often convenient 
to define the component of £ orthogonal to u, 



The results ( |2.26|) and ( |2.27| ) have several interesting features. First, the presence 



k 2 . 



(3.1) 



(3.2) 



The equations of motion of the superfluid are given by 




(3.3) 



13 



In |l6| £ was defined with an opposite sign. 
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together with the constitutive relations 



T» v = (p + P)u»u u + Prf v + f^C + TZ S 
r = qu» - ft" + r dlss (3.4) 

U-£ = /i + fidiss 



where we have chosen to work in an arbitrary 'fluid frame' (see |H| for an explanation 



of this terminology and a fuller introduction to dissipative superfluid dynamics). 

As was the case for the theory of charged fluids which we described in the previous 
section, superfluids also allow for a simple 'canonical' entropy current [TU 



J>i = S u» - - ^4t^ (3-5) 

O canon rri CLISS rp \ / 

where s is the thermodynamical entropy density of our fluid and is related to p and P 
through the Gibbs-Duhem relation 

p + P = sT + fJl q (3.6) 

and 



dP = sdT + qdp+-fd£ 2 (3.7) 
2 



where 



£=V-£^- (3- 



It has been demonstrated in [jilj that the entropy current ( |3.5|) is invariant under field 



redefinitions. It was also shown in [[T(| that the divergence of this entropy current is 
given by 

« = ~d, (|) TZ S ~ (d, (£) - f ) JLs + ^ (/*") (3-9) 
The rest of this section closely follows section |^. In 3A we list the independent 



first order data and second order scalar data, in section |37| we construct the most 
general positive divergence parity conserving entropy current consistent with Lorentz 
invariance. We find that up to a certain ambiguous term which is physically trivial, 
the entropy current agrees with its canonical form ( |3.5|) and therefore the analysis of 
m follows. 



3.1 Onshell inequivalent First order independent data 

In the case of superfluid dynamics, the SO (3, 1) tangent space symmetry at any point 
is generically broken down to SO (2) by the nonzero velocity fields and In the 



- 20 - 



special case that and £^ are collinear, 5*0(2) is enhanced to SO (3). This special 
case is physically interesting since it implies that the superfruid component is motionless 
relative to the normal component — once the superfruid velocity is too large superfluidity 
breaks down. We will find it convenient to decompose all first order fluid dynamical 
data into representations of 5*0(2) and treat the collinear limit as a special point in 
parameter space. 

Representations of 50(2) are all one dimensional. We refer to fluid dynamical data 
that is invariant under 50(2) as scalar data. All other fluid data has charge ±m under 
50(2), where m is an integer. There is always as much +m as — m data. We will find 
it useful to group together +1 and —1 charge data into a two column which we refer 
to as vector data; similarly we group +2 and —2 data together into tensor data. 

Now consider a vector V M whose m — 1 and m = — 1 components are (a,b). The 
vector 

•)> = e^u^Vf, = *V M (3.10a) 

is a pseudo vector. Its components are proportional to (a, —b). Thus, when considering 
representations of 50(2), the same data can be packaged into either a vector or into 
a pseudo vector. The story is similar for all nonscalar representations. For instance, 
a traceless symmetric tensor 7^, whose m = 2 and m = —2 components are (a, b) is 
simply related to a traceless symmetric pseudo tensor 

%u = e tMxM u a ^T rv ee *% v (3. 10b) 

with components proportional to (a, —b). Hence all tensor data can be packaged into 
pseudo tensors. 

We now turn to a listing of the one derivative fluid dynamical and field data for 
superfluids. In Table ^|we explicitly list all one derivative data, one derivative equations 
of motion, and then eventually independent one derivative data. The scalar £ used 
in this table is given by (|3.8| ). We do not list pseudo vectors and pseudo tensors 
independently from vectors and tensors as they are isomorphic and contain the same 
data. In Table |] we assign labels to our independent data. In the same table we 
also present a second listing of a basis for independent scalar data which will be more 
convenient at places. In Appendix |A] we demonstrate that both sets of seven scalars 
and the seven vectors listed are independent data, i.e. that we can solve for all other 
scalars and all other vectors in terms of the chosen basis. 

As can be seen from Tables |3| and f|, after imposing the equations of motion we 
have six first order scalars and one first order pseudo scalar built out of fluid data, 
one first order scalar and one first order pseudo scalar built out of background field 
strengths, five first order vectors built out of fluid data, two first order vectors built from 
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Classification 



All data 



Equations of motion 



Independent 
data 



Scalars (set 1) 



u»d» (#) 



Scalars (set 2) 



u^dpT 



d^=cE^B^ 



(£) 



[d^ v -d v ^)=i-E 
u^8 v T^=-E^ 



P» v dpU u 
P^d^T^) 

«/*e„ (#) 



Pseudo scalars 



Vectors 



P» V E V 

pP u d u T 
P^ u Pd p u h 



P^d v ^ 
pP v i a d v u a 



PP v dpT p u =pP v F uP jP 

P^w {dpi v - d v tp)=P a »Ep 



P^uPdpU v 
PP v uPdpi v 
pP v iPd p u v 

pp v ipdpi v 

PP~ V E V 

Pp»F u ^ 

PP v d v ± 



Tensors 



a 



a 



a 



0". 



Table 3: One derivative data for superfmids. The first column lists all quantities formed 
from the action of a single derivative on fluid and background fields. The second column lists 
all one derivative equations of motion. The last columns lists a choice of independent data. 



The tensors Op\ v and <j" are defined in ( 3.12 )-( 3.11 ). We also used e 



pvafi 



background fields and two independent tensors. The first tensor is simply the usual 
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i 


1 


2 


3 


4 


5 


6 


7 


s? 




P^d^TQ 




eo, (|) 


ed, (*) 




E4 


sf 




P^d^TQ 




u"0„ (J) 


U p dp (#) 




E.tL 


v? 


P» v u p d p u v 


P^u p d p i v 






P pv d v ± 


P pv d v Ep 


p» v F vfi e 



Table 4: Labels for the two sets of independent one derivative scalars and one set of inde- 
pendent vectors. 



shear tensor a^ v projected orthogonal to the plane formed by the two fluid velocities. 

^ = p, a pup ^ _ l^p 75(7 7^ . (3 . n) 

The second tensor c£ is defined by 

oy„. = l -p^p^ (d a ^ + df,£ a - P^P^d^s) . (3.12) 



The counting of data in the absence of background fields agrees with [|10] . 



As was the case for normal fluids, the divergence of the first order superfluid entropy 
current is a sum over quadratic one derivative terms and two derivative pieces of data. 
In order to assist the analysis of the positivity of the divergence of the entropy current 
we list all the scalar two derivative data, the two derivative equations of motion and a 
basis for onshell independent two derivative scalars in Table |5]. Note that we have nine 
independent pieces of two derivative fluid dynamical data (as reported in [IIJ ) together 
with four additional pieces of two derivative data from background field strengths. In 
Appendix [A] demonstrate that the scalars listed in the last column of table [5] form a 
basis of onshell independent scalars. 

3.2 Constructing the entropy current 

With the independent data at hand we proceed with our analysis. The most general 
entropy current allowed by symmetries takes the form 

7 7 7 

j"s = + «" E ^ + ? E «W + E v W> ( 3 - 13 ) 

i=i i=i i=i 

where the coefficients s", s^, and Uj are, at the moment, arbitrary functions of ^, T, 
and £ 2 . Note that we have chosen to expand the terms proportional to u p in the basis 
5" while terms proportional to £ M are expanded in the basis S\. This choice will prove 
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P^u p d p dau u 

H P 1 u 


P^^dpdaUa 

"5 fJ fj, fj, 


P^uPdpdJT^) 


PP u ^dpdp(TCu) 




p»P>'u f 'd a B.,u v 

^ P P V 


s ^p^p \ T J 


uPf p d n d„ (£) 


uP£Pd„d,,T 


All data 


ppua a / p \ 

± OpOy \ 7ji ) 


ppv a fx ( £ \ 
j OpO v 1 y 1 


pp v dpd v T 


u p u pdpdp (|) 






cwfy (£) 


eii p d p dpT 






u^ufdpd^T 




u p^uPdpdpu u 


^PPP-dpd^Uv 








PP y dpE y 


eP^ipU X dpB a 




u^dp (dpJt 1 ) = u 




tfdp (dpJP) = tfdp (cEPB^) 




u 0/3 \S,p,O v l >~ ) - 


- u Op {q^r p U J ) 


{tpd u TP») = 


S% (eFp U r) 


Equations 


vPdf, (UpdyT^) 




(upd v TP») = 


t% {-E^JP) 


of motion 


U^Ob (^U U (dutv 




£% {d£ v - 


-d u £.p))=^dp(eEp) 




d a ( P a ^u v (d u t v 


- d v £ u j) =d a ( P^Ea) 


dp (prdpTfy = 


dp (P^F vj3 J^ 




a ( paptv (X t 


fl C U - /) f pap ?? | 
°vKp) ) — °a\! FpuZ J 






Independent 
data 


P^uPd p dpU u 

^i v upd p dpu v 


P^fpdpdpUp 

u^dpdp (|) 


pp" u pdpdp(T^) 
u^Pdpdp (£) 


PP^Pd p dp(T^) 
uP^dpdpT 




P» v dpd v (£) 




ea u 9pE u 


P^dpE u 




pi™ {dpF vf3 ) ^ 









Table 5: Two derivative scalar data. The first row gives all two derivative scalar data, the 
second row lists all the equations of motion. The third row represents a particular choice of 
independent second order data. 



algebraically convenient below. In total we start with twenty one free parameters in 
the entropy current. 

The two derivative terms in the divergence of the entropy current fl3.13| ) are given 

by 

<Vs = 0? + «i) PP»uPd p d,u u + (s a 2 + v 2 ) PP v uPd p d,{T^) + (s\ + v 3 ) PP^Pd p d,u u 
+ (4 + v 4 ) PP"ed P dp(T^) + (s a 3 + s b 3 ) eCu p d p dpU u 

+ {4 + 4) ^ p dpd, (|) + (4 + 4) u^ p dpdp (|) 

+ (s« + si) uP^dpSpT + v 5 P pv dpd v (I) + s? u P CdpE u 
+ 4 eCdpE u + v 6 P pv dpE v + v 7 (dpF u ?) + . . . . 

(3.14) 

Following the algorithm of the previous section, we first set the coefficient of each of 
the thirteen independent two derivative terms listed in Table |5], which appear in the 
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divergence of the entropy current, to zero. The vanishing of the nine fluid dynamical 
two derivative terms yields the following nine relations between the s a 's s 6 's and u's 

Vl = si v 2 = -s% v 3 = -s\ v A = -s b 2 s b 3 = -sg (3.15) 

S\ = ~8% 8\ = -8% s\ = -S% V 5 = 0. 

The vanishing of the four electromagnetic field related two derivative scalars yields the 
additional four relations 

s* = s b 7 = v 6 = v 7 = 0. (3.16) 

Apart from the two derivative fluid dynamical and background electromagnetic field 
data, there are four nontrivial curvature invariants one can form out of the contractions 
of u^, £ M and g^ u with the the Reimann tensor R a pnu- 14: After plugging in the constraints 
in (|3.15|) and ( 3.16 ) into the expression for the entropy current ( 3.13 ) we find 



+ (Ts a 2 + si) Cu X Rx^ + S a 3 U a ^U^ 6 Ra^6 + ■■■■ 

Each of the terms in ([3.17|) is of indefinite sign. Thus, the coefficients of these four 
terms must vanish. This implies 

sj = 0, s\ = -Ts a 2 , s b 2 = 0, s 3 = 0. (3.18) 

To summarize, by setting the two derivative and curvature terms that appear in 
the divergence of the entropy current to zero we have eliminated 9 + 4 + 4= 17 of the 
original 21 coefficients and are left with an entropy current with four undetermined 
coefficients, 



(3.19) 



The entropy current (|3.19|) can be rewritten in a simpler form by introducing the 
antisymmetric tensor 

Qiiv = T(^u u - e,^) , (3.20) 
introducing a unified notation for the three thermodynamical scalar fields 

^ = ^,4,7'5> / = 1.2.3 



14 We omit the curvature scalar R in this listing since it is a pure gravitational term and therefore 
never appears in the divergence of fluid dynamical entropy current. 
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and also redefining our coefficient functions 

Tfi Ti _ 2 

C -S 2 Ci-S S - ^ 2 _ ^ 2 C 2 - S 4 + ^ 2 _£ 2 C 3 ~ S 6 ~ 7f- 



Then ( |3.19 ) takes the form 

3 

4 = JL^ + coW + J2 CiQTdJli . (3.21) 



i=i 



The divergence of the entropy current ( 3.21|) is given by 

+ (d^c ) (dfjTii) d v Q"> + ad^QTd^i + {d^ Ci ) [dfjTij) QT (d„E<) 



(3.22) 



The first line of the right hand side of ( p.22|) corresponds to the divergence of the 
canonical entropy current and the second line corresponds to the divergence of the 
new terms. The right hand side of (|3.22| ) can be written as a sum of three classes of 



quadratic forms: a quadratic form in one derivative scalar data, a quadratic form in 
one derivative vector data and a quadratic form in one derivative tensor data. Let us 
first focus our attention on the quadratic form in the vector data. All vector terms that 
appear on the right hand side of ( |3.22j ) are linear combination of the six independent 
vector pieces of data, 

P a TV^, ^dvQT, P^VI 1 , P au (d v Ili) (z = l,...,3). 

Independent of the structure of the (as yet undetermined) dissipative terms T^ ss and 
J$ iss , the quadratic form in one derivative vectors does not contain squares of c^Hj, and 
P a p d v Qj iV . 15 It is an easily verified fact that a quadratic form that does not contain the 
square of any given variable is positive if an only if it is independent of that variable. 
Thus, positivity of the divergence of the entropy current requires that the right hand 
side of ( 3.22 ) does not contain any term proportional to P^d^i and P a d v Q^ u . 



The term in fl3.22j ) proportional to P"9 a Ej x d v Q^ u appears with a coefficient 
(q + (9^ c ). Setting this coefficient to zero yields an expression for the q's in terms of 
c 

Ci = -dz i c Q . (3.23) 

Inserting ( |3.23| ) into (|3.21|) we find 



J"s = + d„ fafiT) . (3.24) 



15 This follows because these terms do not appear explicitly on the right hand side of the divergence 
of the canonical entropy current. 
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We have ended up with a one parameter set of entropy currents ( 3.24j) specified by 



the undetermined coefficient Cq. The divergence of these entropy currents is, however, 
independent of c . In fact, due to the antisymmetry properties of Q^ u 

<V£ = <V£ canon • (3-25) 

The parameter cq is essentially trivial and is related to a pullback ambiguity as we 
now explain. It was pointed out in |2B[ that the following set of operations maps one 



positive divergence entropy current J$ to another 

1. Dualize J M to a three- form. 

2. Shift this three-form by its Lie derivative with respect to any vector field 

3. Dualize the resultant form back to a current. 

The end result of this operation is a shift in the entropy current given by (see eq. 6.6 
in[g§) 

5r s = V„( - VJ£) + VVvJZ . (3.26) 

In the current setup we are interested in first order corrections to the entropy 
current. The right hand side of (|3.26|) has an explicit derivative. Therefore, the entropy 
current on the right hand side should be replaced by the perfect fluid entropy current 
Jg = su^. This implies that the second term on the right hand side of ( |3.26| ) is zero 
(recall that the perfect fluid entropy current is divergence free). Moreover V must be 
a derivative free vector field. 

In ordinary (non superfluid) fluid dynamics there is a unique vector at the zero 
derivative order — the fluid velocity u^. Since J M oc w M then V 1 oc implies that 
the first term on the right hand side of ( |3.26| ) also vanishes, and so ( |3.26| ) leads to no 



ambiguity in the entropy current at the first derivative order. 

In superfluid dynamics there exist two zero derivative vectors, and Conse- 
quently ( |3.26|) can be used to generate a shift in the current proportional to d v (c Q Ufl ) 
. We conclude that the freedom to add the total derivative term d u (c Q Q u ^) is precisely 



the 'pullback ambiguity' freeedom described in [26 



Going back to (|3.25|), it follows that the constraints on the dissipative terms T t 



i- 1 " 

diss i 



J diss an d Hdiss from demanding the positivity of the divergence of Js are identical to 
the constraints from the positivity of the divergence of the canonical entropy current 
Ja . Such an analysis has been carried out in [(TO] and we will be carried out in more 

O canon ** V 



generality in the next section. We will not repeat it here. 
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4. Superfluid dynamics without parity invariance 



The main goal of this work is to determine the most general expression for the consti- 
tutive relations for a Lorentz invariant but parity breaking superfluid. We tackle this 
problem using the same algorithm which has been spelt out in the previous sections: 
we determine the most general entropy current which is compatible with the symme- 
tries. The restrictions on the constitutive relations are then determined by demanding 
positivity of the entropy current. 

4.1 Onshell inequivalent data 

As we've discussed in section pO, the listing of one derivative pseudo vectors and 



pseudo-tensors of SO (2) can be constructed from the tensors of SO (2) using the e 
tensor. If V M and T^ v are vectors and tensors of SO (2) then we define the associated 
pseudo vectors V M and pseudo-tensors T^ v through: 

We introduce the following terminology: is said to be the star of V M , and similarly 
*Tnv is said to be the star of 7^,. 

Note that is an SO (2) pseudovector for any vector since the * operation 
projects onto the subspace orthogonal to both ii M and where ^ was defined in 

= P^tv ■ 

In Table |] we list the independent one derivative pseudo-scalars, vectors and tensors 
which we will use in this section. It proves convenient to choose our on-shell independent 
scalars S? to differ from both the choices made in section |[ The linear independence of 
this choice of scalars is demonstrated in Appendix |C[ Our choice of on-shell independent 
vectors is also different from the one is made in Table (| of section |3|. We have used 
the superscript c for the pseudo vectors and vectors to flag this difference. We have 
demonstrated the linear independence of this choice of vectors in Appendix |C|. 

In table [7] we list all two derivative pseudo scalars. The format we use is identical 
to that of previous sections. We list all possible two derivative pseudo scalars that 
can be constructed out of fluid and field data in the left column. The equations of 
motion appear in the second column. For technical reasons we find it convenient to 
formally treat E^ and as independent fields in the full list of two derivative data 
(the left column in table f7|) and to impose the Bianchi identity (which relates some of 
the derivatives of these fields to others) at the same stage that we impose the equations 
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Table 6: Parity odd independent one derivative data which we use in this section. The * 
operation that takes us from vectors V and tensors T to pseudovectors V and pseudotensors, 
T is given in ( [4.1|) . We list a basis for pseudovectors, and a new basis for scalars and vectors 
which we will work with in this section. Note that Vf = *Vf . 

of motion. In other words, in our listings below we will pretend that and 
are independent fields constrained by "equations of motion" which are given by the 
Bianchi identities. In the last column we list a basis of independent two derivative 
pseudo scalars. In Appendix |C| we have demonstrated that the two derivative pseudo 
scalars presented in Table [F] are independent. 



All data 


equations of motion 


Two derivative data 


P^VuB^ 
e X(7 V x E a 
e Xa (u-V)V x u a 
e ACT (C-V)V A M ff 
e ACT (u-V)V A Cr 
e ACT (C-V)V A Cx 


e ACT (CV)C ACT = 
e Xa (V e T° - F ae f) = 
e x °u e V x C 6(T = 
e^ x °^V u F X(T = 
e^ x °u^V u F Xa = 


e Xa V x E a 
e x °((-V)V x u a 



Table 7: Parity odd two derivative scalar data. The first column lists all possible independent 
two derivative scalars with and E^ 1 treated as independent, the second column lists the 
equations of motion (first three rows) and Bianchi identities (last two rows) and the last 
column the independent data. We have used the shorthand notation e A<T = e^^u^Qy and 

4.2 Constraints from vanishing of two derivative and curvature terms 

We denote the total entropy current by 

J"s = JL^ + JL ew + 9AcoQn, (4-2) 
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where Js ncw is a parity odd contribution to the entropy current and d v (co<2^) is a 
residual ambiguity in the definition of the entropy current in the parity preserving 
sector described in section |3|. The most general contribution to Jg nevr is given by 

4^ = E +u "il ~ s & + ^ E ^ • ( 4 - 3 ) 

i=l i=l i=l 

This entropy current has eleven new coefficients, each of which is an arbitrary function 
of T, fi and £. In what follows we will swap the variables /i, and £ with 

U C 2 

and X=7f2- ( 44 ) 

As in previous sections, the eleven parameter set of currents is significantly constrained 
by the requirement that two derivative and curvature terms in its divergence vanish. 
It will prove convenient to rewrite the entropy current in the form 



5 

y e^ p ^ v F pcT + T 2 a 9 cu» + Ta 10 B» 

+ aiVf + a 2 V 2 CM + C [a 3 (u-0 + a A (B-()} (4.5) 



The coefficient functions oi in ( |4.9| ) are linearly related to the coefficient functions in 
( |4.3| ). Although we will never need these relations, for the sake of completeness we 
present them in Appendix 

The first term in (^4.9|) is the divergent-less combination 

V Act V, faTuxC) (4-6) 

with a\ an undetermined function of T, v and x- This expression contributes to the 
entropy current but does not contribute to its divergence. We've observed a similar 
expression in the parity even sector of the entropy current in equation ( |3.24| ). As was 



the case in the previous section, this term is physically irrelevant (it does not contribute 
in any way to the constraints on dissipative terms) and can be ignored in what follows. 

We now proceed to compute the correction to the entropy current ( t4.5|) . It is not 
difficult to verify that the two derivative and curvature contribution to the divergence 
of the entropy current is given by 



V*JLo = «i [^ A XC] V X E„ + (a2+ 2 C " 3) [e^u^u] (CV)V AMc 



D Rp6a\ 

tt P cr\6 — 



(4.7) 



+ 
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where the dots denote the derivative squared terms. The expression on the right hand 
side of ( |4.7p is a sum over a set of linearly independent two derivative data and curvature 
terms (see table |7|). Positivity of the divergence of the entropy current requires that 
the coefficients of the two derivative data and curvature terms vanish. Thus, 

ot\ = «2 = «3 = 0:4 = 0. (4.8) 

Had we not demanded the vanishing of curvature terms above, we would have erro- 
neously concluded that a 2 is non zero. Inserting (fO|) into (|4~5l ) we conclude that the 
most general parity odd addition to the entropy current whose divergence has no two 
derivative data or curvature terms is given by 



2 



c (J, 

5 



+ -rf e^^ u F pa + T 2 a^ + Ta w B» . (4.9) 

Expression (|4.9| ) is the final result of this subsection. 

Note that positivity of the divergence of the entropy current implies a correction 
to the entropy density of the form 

-u^ s = s + a 8 S 2 + atTSt + ^ (c TC) - c T^u-du^ . (4.10) 



Equation ( |4 . 1 U| ) implies that in a parity violating system the entropy density, as defined 



via the Gibbs Duhem relation (i.e., the log of the partition function in equilibrium), 
would receive gradient corrections if Co, <J\ or a% are non zero. 

4.3 Constraints from positivity of the quadratic form 

The divergence of the sum of the canonical entropy current and the the new piece ( (4.9Q 
has no two derivative or curvature terms, and so is a quadratic form in one derivative 
data. In this subsection we will determine the corresponding quadratic form. The 
requirement that this quadratic form is positive will enable us to deduce the restrictions 
on the oVs in Q4.9|) and on the parameters describing the equations of motion of the 
theory. 

From ( |3.9|) and taking into account the anomaly, we find that 

« = -% (|) TZ S +{f~9, (f )) JLs + ^ (/*") + " cvE-B . 

(4.11) 
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On general grounds the divergence of the non canonical part of the entropy current, 
together with the contribution from the anomaly, takes the form: 16 



7 7 



duJZ 

H 1 O new 



(4.12) 



Recall that V Cfl = so that B^ is an antisymmetric 7x7 matrix. 

The expressions for B\, and B\, are straightforward but rather tedious to compute. 
We provide some details relevant to the computation in Appendix 0. We will not 
completely list these matrices here, but will unveil relevant aspects of their explicit 
form as we go along. 

To complete our analysis, we need to reexpress the divergence of the canonical 
entropy current ( |3.9|) in terms of on-shell independent data. We carry out this compu- 



tation in a frame invariant manner, similar to that of section 2.1 



We find it convenient to decompose the explicit fluid first derivative data that 
appears in (|3.9|), namely V M (%) and E^/T — d^v into 5*0(2) invariant scalars, vectors, 
and tensors as follows 

1 



+ 



u-dT 1 



1 

T 



U-VUr 



C 2 T 

d a T 



T 



£ a d a T 
T 



2 1 



(4.13) 



( 2 T 



Ey a 



C-Vi 

Viv + UpU-du + C^-— 



We will find it convenient to organize the 7 scalars and three vectors above into two 
columns and v^, 



udT 

—u-Vv 

_(Vi 

9m (/V) , 
T / 



V 



-V? 

v 1 v 



(4.14) 



/ 



16 This follows from the observation that the final result is a pseudo scalar. The absence of a tensor 
times pseudo tensor term in this expression is not completely obvious but may be checked. 
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The pieces of scalar and vector data in ( f4.14|) multiply particular scalar and vector 
components of T^ ss and J^ iss - It is convenient to classify the scalar, vector and tensor 
combinations as 



Sl 

s 3 
s 5 

t 



rplif p 
± diss r V-v 

^ ' J diss 
fJ'diss 

W ^- L diss 1 a 

P U J a 

pap firp 

r a r v A dissaf3 



C S 2 — (-Tdiss'C 

s 4 = u-T diss -( 

^6 C ' J diss 

V2 = C,T d ZP\ 



(4.15) 



2 



diss,a/3 j 



and to group these into row vectors 



s — (s x s 2 s 3 s 4 s 5 s 6 s 7 ) 
v = (vi v 2 v 3 ) . 



(4.16) 



These definitions permit a very simple formal expression for the divergence of the 
canonical entropy current: 



d J M 

A* S canon 



-s s d 



1 

T 



(4.17) 



In presenting ( [4.1 7|) we have used matrix notation: the expressions on the right hand 
side of ( |4.17|) are each the product of a row and a column. 

Equation ( |4. 17| ) is not our final expression for the divergence of the canonical 
entropy current for an important reason; the entries in the columns and are not 
independent on-shell. In fact, on-shell, each of the 7 is a linear combination of only 
four independent scalar terms. Similarly, on-shell, each of the three entries in is a 
linear sum over two on-shell vectors. Specifically 



A V V 



(4.18) 



where 



S 



Si 



V 



(4.19) 
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and 



A s 



Rs 

2g n T x 
Rs 



R_ 



1_ 





B 3 A 3 


B 2 A 2 


Bi Ai 


3T 2Tx 


3T 2Tx 


3T 2Tx 


B 3 i A 3 
3T T T X 


B 2 1 A 2 

3T T X 


Bi , Ai 
3T ^ T X 


1 

T 2 










K 2 




T 


T 


T 





-1 














{p+P)K 3 


(p+P)K 2 


(p+P)K 1 


T 


T 


T 



A v 



f-R 




V 



2 

1 



(4.20) 



with 



R 



(4.21) 



p + P 

Note that S is a column composed of 4 of the 7 on-shell independent scalars listed in 
table ^ while V is a column composed of 2 of the 6 on-shell independent vectors listed 
in the same table. Note also that A s is a 7 x 4 matrix while A" is a 3 x 2 matrix. The 
computation of the entries of the matrix A s and is outlined in Appendix |C| where we 
also provide expressions for Aj, Bi and i^j. It follows that when reexpressed in terms of 
on-shell independent data, the divergence of the canonical part of the entropy current 
takes the form 



d J M 

A* S canon 



-sA s S - v^V" - ^ V7? 



•/if 



(4.22) 



Equation ( |4.22| ) depends only on those 4 linear combinations of the 7 scalars Sj 
that appear in the 4 columns of sA s . Similarly ( f4.22| ) depends on only the 2 linear 
combinations of the three vectors Vj that appear in the 2 columns of \A V . These 
numbers could have been anticipated on general grounds. As in the case of ordinary 
fluid dynamics (see section §) the equations of superffuid dynamics suffer from a field 
redefinition ambiguity. The field redefinitions in question are T and /x, and may 
be decomposed into 3 SO (2) scalars and one SO '(2) vector. All of the seven scalars Sj 
described above transform under the scalar field redefinitions. As these redefinitions 
depend on three parameters, 7 — 3 = 4 independent linear combinations of the scalars 
are field redefinition invariant. Similarly we expect to find 3 — 1 = 2 field redefinition 
independent linear combinations of the three vectors Vj. Since the canonical entropy 
current is field redefinition invariant, the expression for its divergence is necessarily 
also frame invariant. It follows that SiA^ and ^iA v - must be linear combinations of the 
four frame invariant scalars and the two frame invariant vectors respectively. We have 
explicitly checked (see Appendix 0) that the the four linear combinations sA s and the 
2 linear combinations \A V are indeed field redefinition invariant. 
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Putting together fl4.12j) and (|4.22|) , positivity of the entropy current implies that 



-^V7T > (4.23a) 

-vf^.v; + vr-^.v;>o (4.23b) 

- Sl A^ + ^B^>0. (4.23c) 

In the rest of this section we will deduce the constraints imposed by ( |4.23| ) on the 
one derivative constitutive relations for the tensor t^, the frame invariant vectors v^A v 
and the frame invariant scalar s sA s . 

4.3.1 Positivity in the Tensor Sector 

The most general one derivative constitutive relation for the frame invariant tensor t 
(on symmetry grounds) is given by 

^•[iv Py Tdissa/3 O^"^ ^ Pdiss,otf3 ^]{T% p,v ^jiTifjiV (4.24) 

where T\ = <J U , % = and % = *Ti and the 77's and rfs are four possible coefficients. 
Positivity of the entropy current requires that 

r)i > 771 G R r]2 = V2 = 0. (4.25) 

In other words, the tensor part of T£" ss has no term proportional to T2 and to *72- The 
latter result could have been anticipated. The divergence of the entropy current does 
not have a piece proportional to T2 ■ Therefore terms in the divergence of the entropy 
current which are linear in 7~2 must vanish. On the other hand does, in general, 
have a term proportional to 7i; the coefficient of this term is the shear viscosity of the 
fluid and is constrained to be positive. Positivity also allows a term proportional to 
*T\. The coefficient of this term, fji, is a new parity odd coefficient (it is of course 
an arbitrary function of the thermodynamical scalars). Note that fji drops out of the 
formula for entropy production. This is because the contraction of any tensor with its 
own star vanishes. Consequently, fji is non dissipative; its coefficient is unconstrained 
by the requirement of positivity of entropy production. 

4.3.2 Positivity in the Vector Sector 

The most general constitutive relations allowed by symmetries for the two frame in- 
variant vectors takes the form 

v?4£ = -V^ - Vf^ . (4.26) 
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The matrix of possible transport coefficients k is a 7 x 2 matrix. The same is true for 
k. The first index runs over a basis of on-shell inequivalent vectors, while the second 
index runs over a basis of field redefinition invariant combinations of vectors in T^ ss 

and J diss- 

A useful observation is that the quadratic form in the vector sector contains no 
terms proportional to the square of Vf for i = 3 ... 7. The first expression in ( (4.23b| ) 
has an explicit factor of either Vf or Vf. Since is antisymmetric, the second ex- 
pression in ( |4.23b|) has no terms proportional to the square of any of the on-shell 
independent vectors. Positivity thus demands that the quadratic form in the vector 
sector be completely independent of Vf for i = 3 . . . 7. This implies that 



Bl = 


i — 3, . 


• ,r, 


j 


= 3,. ..,7 


(4.27a) 


%ij = —B"j 


2 = 3,. 


-,7, 


j 


= 1,2 


(4.27b) 


f^ij = 


i = 3, . 


• ,r, 


j 


= 1,2 


(4.27c) 



We now proceed to investigate the implications of ( |4.27 ). 

Let us first focus on (|4.27a|) . A priori, this equality gives us a set of 5 x 4/2 = 10 
partial differential equations that constrain the six independent cr^ coefficients in ( [4.9|) . 
One of these equations is a tautology. The remaining nine equations are as follows 

d(aj - as) 1 da 3 dcr 5 1 da 3 

B ^~ ^ &T T^7"° 5 35-o^ — --— _o 

w - n =v ~ ag ) 9(75 - n r» - n ^ ^(^io - ^s) 03 _ n 



™ n ^ g((Ti -^a 8 ) , n ri , 9/t „ n 

# 57 = U =>- h cr 5 = # 46 = =^ — 2<Jio + 2z/cx 4 = 

ax 



BJ 7 = ()=*• ^-^) +(j4 _ Jg) = c „ (4 28) 

Equations (|4.28 ) are far from independent. They admit the following two parameter 
set of solutions 

,0(010 - i/0- 8 ) 



03 = 



dT 

d(a w - va H ) 



04 = 08 + cv - 

a, \ 6V (4.29) 

°" 5 = 5 

2 , 

09 = s 9 - 7TCZ^ + 2l/(cTio - Z/0s) 
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In ( 4.29Q we have chosen as and <7iq as our free parameters. The term sg in the last 



line of ( 4.29|) is an integration constant. It turns out that the requirement of CPT 
invariance forces 

s 9 = 0. 

The argument for this conclusion is very similar to that given below ( p.27| ). While 
we present all formulae below at nonzero Sg, the reader should keep in mind that sg 
actually vanishes in any superfluid that enjoys invariance under CPT invariance, i.e. 
any superfluid that arises from a quantum field theory (and so presumably for any 
superfluid in the real world, or obtained via the AdS/CFT correspondence). 

Equation ( 4.29| ) is the most general solution to ( }4.28| ). As an aside we note that 



(|4.29|) also implies that B\- = for i, j = 1, 2. We will use this information shortly. 



We now turn to the implications of ( f4.27b| ). In this case, the coefficients k whose 
second index lies between 3 and 7 are completely determined in terms of By, which is, 
in turn, a function of the two free parameters a 8 and a w . We find that B i2 = for 
i = 3 . . .7. Thus, 

k i2 = i = 3, . . . , 7. 
However kn are in general nonzero. They are given by 



K 31 — 


-B v 31 = 


—RTcr 3 — T8tcT8 


k 41 = 


--B41 = 


-RT 2 a A - Td„a$ 


«51 = 


-Bt % = 


-RT 2 a 5 - Td x a 8 


kei = 


-#61 - 


-2RT 3 a 9 + 2T 2 a 10 


k 71 = 


TDV 

-B n - 


-RT 2 a 10 + 2Tcr$ + cTv 



(4.30) 



To be clear we reiterate (|4.27cj ) 



% = % = 3,..., 7, j = 1,2. (4.31) 

Once we have implemented all these conditions, the quadratic form in the vector 
sector takes the form Vf-^Vj + V^-k^Vj where i,j = 1,2. Positivity of the entropy 
current implies 

E E ( v <' + v?-*w) > 0. (4.32) 

i=l,2 j=l,2 

The Onsager relations (see [0) imply that K12 = K21 and k±2 = —k 2 i. We are left with 



6 independent transport coefficients Ku, k 2 2, K12 = ^21, ^22 and k \ 2 = —R21, each 
of which is an arbitrary function of /x, T and ^, but is constrained by the inequality 
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that the matrix 



K 



(4.33) 



Ku (k 12 Ki2 
K 12 (k 12 «22 
\~C«12 «12 K 22 y 

is a positive semi-definite matrix. Owing to the high degree of symmetry of /C, this 
condition reads 

KH«22 > «? 2 + C 2 ^i2 • ( 4 -34) 
Note that the 3 k transport coefficients occur even in parity preserving superffuids, 
while the 3 k coefficients violate parity and are new. Note also that two of the three 
new parity violating transport coefficients, kn and k 22 , do not appear in the expression 
for entropy production, and so are non dissipative. However the third new coefficient 
k\2 causes dissipation. 

4.3.3 Positivity in the Scalar Sector 

The most general expansion of the constitutive relations in the scalar sector takes the 
form 

7 2 

s ^'i = - E - E > ( 4 - 35 ) 

i=l i=l 

implying the inequality 

7 4 2 4 2 7 

E szfasj + E E + E E * w ^ • ( 4 - 36 ) 

i=l j=l i=l J=l i=l jf=l 

Terms in ( |4.36|) involving products of scalars and pseudo-scalars must vanish since 
they can not be positive definite for an arbitrary flow. Thus, focusing on the parity 
odd sector we must impose 

2 4 2 7 

E ^ s » s i + E E s c i BS n s i = • ( 4 - 37 ) 

i=l j = l i=l j = l 

Only the second piece on the right hand side of ( |4.37| ) involves Sj with j — 4 . . . 7. We 
immediately deduce the equation 

fl?. = =1,2, j=4,...,7. (4.38) 

Equations ( |4.38| ) are automatically satisfied once we impose the solution (|4.29|) ; we have 
no further restrictions on the free parameters a 8 and a w . For j = 1, ... ,4, equation 
( )P7D implies 

h = -Bit ■ (4.39) 
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It follows that there are no undetermined parity odd transport coefficients in the scalar 
sector; all parity odd contributions to frame invariant scalar combinations of T^ ss , J^ iss 
and Hdiss are completely determined in terms of as and a\o by ( f4.39| ) and the explicit 
listing 

B s = ( ~ " 2fT3 " 2T2R ^ - 2TfJ 4 " 2T 2 K 2 a 9 -2Ta 5 - 2K x T 2 a<\ 

(4.40) 

Note that B\ x = —Rqi/( 2 and = —K, n /( 2 with K a defined in ( [4.30| ). 



The remaining undetermined constitutive relations in the scalar sector are associ- 
ated with the parity even coefficients. These are parameterized by fyj where i and j 
both range between 1 and 4. The Onsager relations imply that the antisymmetric part 
of flij is zero. So we have a total of 10 parity even dissipative coefficients in the scalar 
sector. These 10 coefficients parameterize a symmetric matrix that is constrained to 
be positive 

4 4 

^^S<7V>J>0. (4.41) 

i=l j=l 

5. Summary of results and special limits 

In sections [3] and ^ we have built up a complete theory of relativistic superfluid hy- 
drodynamics that may or may not enjoy a parity invariance. Our work generalizes the 
recent work of |10| (see also |L6]) for parity conserving superfluids, and the work of |20 



for charged gauge-theory fluids with triangle anomalies (see also [^T], [2~2]j). 

The computations presented in sections |3| and |] are lengthy and algebraically rather 
intensive. In this section we present the final results of our calculations in a manner 
that is self contained and makes no reference to the method of derivation. 

5.1 The problem addressed 

As discussed in section ^, the energy momentum tensor T^ u and charged current of 
an s wave superfluid are given by 

= ( P + p)u»u u + p^ v + fee + tz s 

J" = q n u» - fe + JLs (5-1) 

= /i + Vdiss 

where £^ is related to the background gauge potential and Goldstone Boson ip 
through 

^ = -d^ + A,. (5.2) 



- 39 - 



(ip is the phase of the charged scalar condensate.) We often use the projected variable 



^ = ^ + uV)^ = P^^. (5.3) 
in what follows. The equations of motion of superfluidity are given by 

V = cEpB" (5.4) 

where we have allowed for the presence of triangle anomalies via the non conservation 
of the charged current. In writing ( |5.1| ) we have not specified a frame — a particular 
definition of the temperature T, chemical potential fj, and four-velocity w M once we 
deviate from thermodynamic equilibrium. The thermodynamic functions p, P and / 
are functions of the variables 

T, v=£, X = ^f- (5.5) 

The equations Q5.4| ) and ( |5.1| ) together specify a complete set of equations for the 9 
fluid dynamical fields m m (x), T(x), h(x), once T^ ss , J% iss and fidiss are specified 

as functions of the fluid dynamical fields and their derivatives. The equations that 
determine T£? s , J^ iss and fidiss in terms of derivatives of fluid dynamical fields are 
called constitutive relations. The main result of this paper is the determination of 
the most general form of the constitutive relations consistent with Lorentz invariance, 
positivity of the divergence of the entropy current and the Onsager relations. 

The most general allowed constitutive relations turn out to depend on 20 unspeci- 
fied functions of T, v and % which we will refer to as free parameters. In this section 
we will present our final results for the constitutive relations in terms of these 20 pa- 
rameters. 

Superfluid flows are also accompanied by an entropy current. The requirement that 
this entropy current be of positive divergence for every fluid flow played a key role in 
our derivation of the allowed form of the constitutive relations. At first order in the 
derivative expansion the entropy current takes the form 

J S = J% +J% . (5.6) 

o canon o new \ / 

where Ja is the so called canonical entropy current (in a fluid frame) 

O canon \ / 

J Scanon ~ iU rr J dlSS rp \° • ' J 
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and Js ncw is a correction proportional to a sum of single derivatives of fluid fields. In 
this section we also present our results for the correction J a to the entropy current. 

1 O new 1 m 

After summarizing our results we proceed to explain the specialization of these 
results to the case of conformal or Weyl invariant superfluids, and certain simplified 
but often used limit where the superfluid velocity is relatively small — in liquid Helium 
once the superfluid velocity is too large superfluidity breaks down. 

5.2 Listing our main result 

The most general correction to the entropy current is given by 



JL, = ^ pa 9 v (<tiTu p C) + d u (coQ v n + a 3 V c / + Ta,V c / + Ta 5 V 5 



c n 



where 



Q/J.U = T(^u u - , (5.9) 
and cr 3 , <7 4 , <7 5 and a 9 are determined in terms of a 8 and a w by 

cr 3 = -T-|^(crio - vcr 8 ) 

d d » (5.10) 

dx 

2 

a 9 = s 9 - -cv 3 + 2is(a 10 - z/cr 8 ) 



In the above relations u and B were defined in (|5.24|) , c is the anomaly coefficient and 



Sg is an integration constant. In any CPT invariant theory 

s 9 = 0. 

The terms in the entropy current proportional to Co and o\ are divergence free and 
so are physically irrelevant (they have no effect on the equations of motion). Ignoring 
these terms, the final result for the entropy current is expressed in terms of as and crio 
which are undetermined functions of T, v and x- As we will see below, as and <Tio enter 
the constitutive relations, and so are free parameters. 17 



17 



In a parity preserving theory all the <7i 's are forced to vanish. As a consequence, the correction to 



the entropy current in a parity even theory is restricted to the second term in (5.8) which is physically 
unimportant. For many practical purposes the entropy current in parity even theories may be taken 
to be uniquely fixed to its canonical form. 
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We now turn to our results for the constitutive relations. We express our results 
in a field redefinition invariant manner. Since the expressions T^ ss , J^ iss and ^iss are 
not separately field redefinition invariant, we first list field redefinition invariant linear 
combinations of these quantities. Let 



where 



Sl 

s 3 
s 5 

t 



rpflV p 
± diss r 

U ' Pdiss ' W 
^ ' J diss 

pdiss 

T m pv 
a ^ ± diss r a 

P\J a 



pap firp 



diss a/3 



_ _ pV v potlirp 

„ r r 1 diss,af. 



S 4 

s 6 



C ' Pdiss ' C 
^ ' Tdiss ' C 
C ' J diss 

(• rplMCt pv 

>M diss a 



We define the row vectors 



if + u^u 



pilV _ pfjLU _|_ 



(5.11) 



(5.12) 



S — (si S2 S3 S4 S5 S6 S7) 
V = (vi V 2 V 3 ) . 



We also define the matrices 



A s 



/ 



Rs 

2g n T x 
Rs 

<lnXT 



_ R 


1_ 





B 3 A 3 


B 2 A 2 


Bi Ai 


\ 


3T 2T\ 


3T 2Tx 


3T 2Tx 


B3 , A3 
3T T T X 


B 2 1 A 2 

3T T X 


Bi , Ai 
3T T X 






1 










K 3 




A'i 




T 


T 


r 







-1 


















(P+P)K 3 


(p+P)K 2 


(p+B)A-i 


J 


T 


T 


T 



A' 



(-R \ 

^ 

V-i 0) 



(5.13) 



(5.14) 



where 



i? 



1 T 



P» v d v v 



P + P 

and the Aj's -Bj's, Cj's and fQ' s defined in appendix f|. 

In terms of ( 5.11| )-( |5.14j ), the frame invariant scalar, vector and tensor combinations 
of T^ ss , J% iss and fx d iss are given by the row vectors 



sA s 



(5.15) 
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i 


1 


2 


3 


4 


5 


6 


7 


i 


P^V 1(T 




P^ a V a T 


P^V a (£) 






7~>IJ.l/ T7I >rv 

FuaC 


s 




u-dT 




























5, 

















Table 8: The basis in which we present our main results. The pseudo vectors and pseudo 
tensors V and f are given through W = e^u^ a Vp and f ^ = e w °>P Up £ a JJ + e u P a Pu p £ a T/f 
Here P> = ^ + it^u" - C M C7C 2 - 



By scalars, vectors and tensors we mean expressions which transform as spin 0, ±1 and 
±2 representations of the 50(2) symmetry that is left invariant by the two vectors 
and £ M at each point in spacetime. 

We have 4 frame invariant scalars, 2 frame invariant vectors and one frame invari- 
ant tensor. The constitutive relations determine these quantities as functions of first 
derivative fluid expressions through the equations 

\r = -vrr - vtr 

2 2 / 7 

v i A l E v ' i; 'j - E v,i; 'i ~ 6 r- ( E V ' /V ' ! 

i=l i=l \ i=3 

4 4 / 2 4 

*^ = - E E ^ -EE ^ 

1=1 jr' = l \i=l 3=1 

with T, T, V, V, 5 and 5 a basis of onshell independent 50(2) invariant tensors, 
vectors and scalars given in table || In equation (|5.16| ) K and are symmetric 
matrices and R12 = — «2i- Those coefficients that occur in the big brackets on the right 
hand side of that equation, namely Ru for i — 3 . . . 7 and for i = 1,2 and j = 1 . . . 4, 
are not free but are determined in terms of a 8 and a w through the equations 

R 3 i = -RTcr 3 - Td T a 8 
K41 = — -RT 2 <t 4 — Td u a 8 

K 51 = -RT 2 a 5 -Td x a 8 (5.17) 
k 6 i = -2RT 3 a 9 + 2TV 10 
«7i = —RT 2 a\Q + 2To"8 + cTz/ 
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-Pi- 



_ 2^0 _ 2fT3 _ 2T 2 K 3 a 9 -2Ta 4 - 2T 2 K 2 a g -2Ta 5 - 2K l T 2 a 9 \ 
T^-T^ + ^f 1 d T a 8 -K 3 Ta 10 d„a 8 - K 2 Ta l0 d x a 10 - K x To XQ ) 

(5.18) 



The remaining coefficients in ( |5.16| ), those that occur outside the big brackets on 
the right hand side of ( |5.16 ), are free parameters. The symmetric 2x2 matrix k, 4x4 



matrix /3, 77 and K i2 = — R 2 \ are constrained to obey 

T] > 

kuk 22 > k 2 12 + C 2 k 2 12 (5.19) 
Pij is symmetric positive semi-definite . 

The remaining parameters 

R u k 22 fj a 8 <7 W (5.20) 

are unconstrained. 

In the introduction we counted twenty free parameters. Let us enumerate them 
explicitly. The parity even parameters consist of a 4 x 4 symmetric matrix 0^ in the 
scalar sector (10 parameters), a 2 x 2 symmetric matrix in the vector sector (3 
parameters) and the shear viscosity 7] in the tensor sector (1 parameter) yielding a 
total of 14 parity even free parameters. These are the only parameters in a theory 
that conserves parity. The parity odd parameters consist of the three parameters k u , 
k 22 , k 21 = —K12, in the vector sector and the one parameter fj in the tensor sector, 
amounting to a total of 4 free parameters. The two undetermined functions a 8 and <7i 
are two additional free parameters. 

Five of the six parity odd free parameters, namely Ru, R 22 , fj, a 8 and aio are 
non dissipative, in the sense that they do not result in entropy production. The one 
derivative fluid data that these parameters multiply can be nonzero even in stationary 
solutions. On the other hand the remaining parity odd parameter R 2 \ = —R\ 2 is 
dissipative, and this is also true of all 14 parity even constitutive parameters. These 
parameters cause entropy production. The requirement that these parameters produce 
rather than destroy entropy constrains them to satisfy ( |5.19| ). 



After completing the work reported in this paper we received the preprint [33] 



which addresses the same question dealt with in this section. While we have not yet 



performed a detailed comparison between our results and those of |33fl , a superficial 
comparison of final results reveals several differences. We list some of these differences 
here. 



1. The author of ]33| finds a unique entropy current, while the entropy current in 



this paper has a two parameter (function) freedom. 
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The author of [33| finds a total of 2 new parity odd dissipative parameters, while 



we find 6 such parameters in this paper. The parameter we refer to as rj seems 
unaccounted for in 



The author of |33| reports that his new parameters are all functions of a single 
thermodynamical variable, while the six new dissipative parameters in this paper 
are all unconstrained functions of all three thermodynamic scalars. 



The author of |!3] finds no mixing between parity even and parity odd sectors, 



while in this paper we find nontrivial mixing between parity even and parity odd 
vector contributions. 

We do not yet understand the reason for the discrepancy between our results and those 



of [33]. We leave an investigation of this issue to future work. 



5.3 Weyl invariant superfluid dynamics 

The equations of superfluidity simplify somewhat in conformal field theories, i.e. theo- 
ries that enjoy invariance under Weyl transformations. The source of the simplification 
is twofold: first the stress tensor of a CFT is always traceless. Second, the dependence 
of all thermodynamic functions on temperature can be deduced by dimensional analysis 
if we take the remaining variables to be the dimensionless quantities v and \. In the 
rest of this subsection we briefly outline the special simplifications in Weyl invariant 
superfluid hydrodynamics. 

Let us first consider the entropy current. The coefficient a 3 in the expression (|5.8|) 
for the entropy current of a Weyl invariant theory must vanish. This follows because V3 
does not transform homogeneously under Weyl transformations. Since as and <7io are 
dimensionless and therefore independent of T, equation ( |5.10|) does not pose an extra 
restriction on these two parameters. 

Let us now turn to the consequences of Weyl invariance on constitutive relations. 
It turns out that (beside implying that all constitutive parameters are independent of 
T) Weyl invariance implies no special simplifications for constitutive relations in the 
tensor sector. In the vector sector we expect that = since, like 03, it multiplies 
V3 which does not transform homogeneously under Weyl transformations. Since 03 
vanishes and <7g is independent of the temperature, this condition is automatic from 
equation (|5.17|) , i.e. it does not pose an extra restriction on as- The remaining terms 



in the vector sector remain unchanged. 

In the scalar sector dimensional analysis implies that 



A 3 = K 3 = B 3 = ~. (5.21) 
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The fact that the scalar S% does not transform homogeneously under Weyl transfor- 
mations implies that both indices of /3 in the third equation in ( |5.16|) run only over 
the indices 1,3,4. Consequently /3 is effectively a 3 x 3 symmetric matrix (with 6 pa- 
rameters) in the Weyl invariant case. More importantly, the tracelessness of the stress 
tensor implies that the number of scalar terms in T^ ss , J^ iss and fi^iss is six rather than 
seven. In other words, the seven terms Sj listed in ( |5.11| ) are restricted by the linear 
relation 

si + s 2 -s 3 = 0. (5.22) 

Since there are still three field redefinitions in the scalar sector, this implies that we have 
three rather than four field redefinition invariant scalar terms which can be constructed 
out of T^ ss , J^ iss and /i diss . Accordingly we find that 

: 2 (5.23) 

#2 = 

so that the row vector SiAfj has three rather than four independent entries. Equation 
( |5~2D follows from ( ggg ), ( ggg ) and a 3 = 0. 

It follows that Weyl invariant superfluid hydrodynamics is governed by 10 consti- 
tutive parameters (when it preserves parity) and 16 constitutive parameters (when it 
violates parity). 

5.4 The C ->■ limit 

In liquid Helium, superfluidity breaks down once the velocity of the superfluid com- 
ponent is too large relative to the velocity of the normal component. The projected 
superfluid velocity = P^^ can be thought of as a control parameter for that limit. 
Once ( = one can locally boost a fluid element to a frame where neither its normal 
component nor its superfluid component are in motion. 

When ( is non zero, the velocity of the normal component w M together with ^ break 
the 50(3, 1) symmetry to 50(2). In the £ — >■ limit (the collinear limit) this symmetry 
is enhanced to 50(3). The number of possible expressions for 50(3) invariant scalars, 
tensors, vectors and their parity odd counterparts is significantly smaller than that of 
the 50(2) symmetric case. As discussed in [|T5|, [TTJ there are two possible scalars, d-u 



and <9 M (f^) (which we can swap with u-dT and u-dv via the equations of motion), one 
vector Vi = —P^dvV + E^/T and two tensors in the parity even sector, o^ v and the 
symmetric traceless projection of d^ u onto the space orthogonal to u^. In the parity 
odd sector we have two pseudo vectors 

^ = l -^ vpa u v d p u a and B»=^ vpa u v F pa . (5.24) 
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Thus, in the collinear limit the transport coefficients of the theory reduce to a 2 x 2 
symmetric matrix ^ (in the parity even scalar sector), a diffusion coefficient k (in the 
parity even vector sector), a shear viscosity rj (in the parity even tensor sector), and 
two parity odd terms k^ and Kb- These enter the constitutive relations through 

d-u mta , „ u-dT T Udi 



o T T dL P nv 7j^- u - T diss-u + u-dvu-J diss + -j^df, (/f) = SiPijSj 



(RP» a T™ s u u + P» a r) = kV? - k B B» - k^ 



P/l T)V rpO/3 1 ryfiv JD rpO.fi 

ol (3 diss ~q ^a/3-l-diss — ~ r l a 



(5.25) 



It must also be true that 



J% new = W + * B B» + ... (5.26) 



where we have written only the odd part of the entropy current, omitting the total 
derivative term in the parity even sector (the second term on the right hand side of 
(|5.8|) ). This omission is represented by the ... in ( |5.26| ). 

As we have explained, we expect ( |5.25| ) and ( |5.26| ) to hold on general grounds. 
We should, therefore, be able to verify these equations — and read off the values of Kg, 
k~ u — by studying the ( — » limit of our general results. We now describe how to take 
this limit. 

Let us first consider the entropy current. The third, fourth and fifth terms in ( |5.8| ) 
each have an explicit factor of ( and so simply vanish in the ( — > limit. (We assume 
that all physical quantities like the entropy current — and therefore all all coefficient 
functions that appear in the entropy current — are analytic functions of (.) The last 
two terms in (|5.8|) are already proportional to and B^ respectively. The sixth term 
in ( |5.8|) is proportional to B^ as £ — » 0, while the first term has a piece proportional to 
u^ 1 and a piece proportional to B^ . Putting it all together we recover the form (|5.26|) 
with 

a u = T 2 (s 9 - ^cz/ 3 + 2u (a w - vo % - -a A ) 

, x n (5-27) 



OB — T 1^(710 - ZAJg - -CTi 

where all functions are evaluated, of course, at ( 2 = 0. 

Let us now turn to the constitutive relations. Notice that — (vA v )i is simply the 
5*0(2) vector part of the SO (3) frame invariant vector presented on the left hand side 
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of the second equation in ( 5.25 ). Note also that at leading order in small ( 



2T(sA s ) 2 2T{sA s ) 2 2T{sA s ) 1 2s 2 - Si 
A 3 A 2 A x X 

Since the constitutive relations ( |5.16| ) equate this combination to an analytic function 
of (, it follows that in the ( — > limit 2s 2 — Si is frame invariant and, in fact, vanishes 
at ( = 0. 18 Using this fact it follows that — (sA s )x is proportional to the scalar part of 
the vector on the left hand side of the second equation in ( |5.25| ). Thus, 

k > 
V > o 

0ij is positive semi-definite /g 28) 

4 

~ Ku] = -k 6 i = -C 2 /?ii = 2RT 3 s 9 - -cRT 3 u 3 - ART 3 v 2 a$ + 2T 2 (2RTv - l)a w 
k b = ~k 62 = ~( 2 (3 2 i = RT 2 cr w - 2Ta 8 - cTv . 

It was necessary that — k 61 = — ( 2 f3u and that k 62 = —( 2 (3 2 \ in order that our results 
group into an SO (3) vector in the ( — > limit. This gives us a mild consistency check 
on our results. 

Let us summarize. The parity odd contributions to the entropy current and con- 
stitutive relations are much simpler in the collinear limit than in the general case. The 
only constitutive relation that receives corrections in this limit is the expression for 
the frame invariant vector listed in the left hand side of the second equation of (|5.25|) 



(equal to the dissipative part of the charge current in a Landau like frame). The cor- 
rections to both this vector and the entropy current are linear combinations of B^ 1 and 
uA These linear combinations are specified in terms of 3 free functions, erg, ct\q and 
<j\ and one integration constant (which vanishes on assuming CPT symmetry). As 4 
coefficient parameters are determined in terms of 3 parameters, the coefficients obey a 
single linear identity. This identity is given by 

1 /i 3 1 2 

-Ou - = -t^c + -s 9 T (5.29) 
Under the eminently reasonably assumption of CPT invariance Sg = and 

-o u -iMj B = -—c (5.30) 



18 This could have been anticipated on general grounds, as this combination is proportional to the 
SO (2) scalar in the SO(3) stress tensor, i.e. it is proportional to (■<?■(. 
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The relation ( |5.30|) is the only real prediction in the SO (3) invariant sector. Note that 
it constrains coefficients in the entropy current alone. In the next section we will use the 
AdS/CFT correspondence to test this relationship (of course Sg = in the holographic 
context). 



The ( — > for Weyl invariant superfluids 

As we have seen above, Weyl invariance removes one of scalar field redefinition invariant 
combinations of one derivative corrections to the constitutive relations (because the 
trace of the stress tensor is set to zero) and also removes one of the scalars in terms of 
which the parity even part of these constitutive relations is expressed (because u ■ dT 
and d-u do not transform homogeneously under Weyl rescalings). In more detail, since 
in a conformal field theory, q n (J^ In u-dv = there is only one permissible 

on-shell scalar. In the ( — > limit the 2x2 matrix described above collapses to a 
lxl matrix. Thus Weyl invariant superfluid dynamics is characterized by 3 constitutive 
parameters in the ( — > limit, one in the tensor sector (the shear viscosity), one in the 
vector sector (conductivity) and one in the scalar sector (a new coefficient). 

Apart from determining the dependence of k^, Kg, and <tb on the temperature 
(using dimensional analysis), Weyl invariance does not impact the discussion presented 
above for parity odd corrections to constitutive relations and the entropy current. In 
particular (|5.25|) - (|5.30|) continue to hold, and do not significantly simplify in a theory 
that enjoys Weyl invariance. 



6. A holographic computation 



The AdS/CFT correspondence JJ5], |36j, |37j can be used to study the superfluid phase 
of appropriate large N gauge theories in the limit of infinite t' Hooft coupling |38|, [39], 
34, [l(| 41]]. We will refer to such phases of gauge theories as holographic superfluids. 



Various features of holographic superfluids have been studied in the literature. These 
include sound modes j|2], [43], [41], [44], |i"6|j , critical superfluid velocities |34], [45[] and 
vortex structure |^6j, 47, [IS]]. Holographic superfluids with a non uniform temperature, 
chemical potential and velocity fields have been recently constructed in [IB, [TO] (see 
also ||). These dynamical superfluids were constructed for a certain specialized bulk 
action which, following the work of [pi, allows for an analytic treatment. 



In this subsection we will demonstrate that the prediction ( p.30| ), that follows from 
our general analysis, is indeed true for superfluids that admit a dual description via the 
AdS/CFT correspondence. Our demonstration uses the technology of the so called fluid 
gravity correspondence 



H |5], H, |7|, m, P, ||, |D but allows for the implementation 
of this map in an abstract manner that does not require us to know the explicit form 
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of the solutions dual even to stationary fluid flows. In this respect our analysis is 
reminiscent of the thermodynamical analysis of Sonner and Withers ||. 
Our starting point for the analysis is the bulk action 



S — Seh + S,, 



'matter 



+ s t 



cs 



(6.1) 



where 



Seh 

^matter 

Scs 



1 
1 

C 

24 



-9 (i? +12) 

-a \ ~V F (\i>\)F mn F mn - V4\ij\)\d m tP - iqA m ^\ 2 - V(\^\) 



_ mnpqr a rp rp 
i/ c ■ rL m- L np ± qr 



(6.2) 



with e 01234 = l/y/=g, F = dA, V F (0) = 1, V^(0) = 1 and V(0) = 0. We will often use 
D m ip = d m ip — iqA m ip. Roman indices run from to 3 and 5. Later we will use Greek 
indices to denote the boundary coordinates \i = 0, . . . , 3 and i,j = 1, . . . , 3 will denote 
the spatial coordinates along the boundary. We have chosen the action ( |6.1| ) for several 
reasons. As has been discussed in detail in |3S], 311, to construct a holographic 

superfluid one needs to spontaneously break a gauge symmetry in an asymptotically 
AdS geometry. The minimal required fields for such a construction involve a gauge 
field Ap, a charged scalar ip and a metric g^. The parity conserving part of the action 



(|6.1|) , i.e., Seh + S ma tter, is the most general action one can write down which involves 
the fields described above. The parity violating sector consists of a Chern-Simons term, 
Scs, which naturally appears in consistent truncations of type IIB supergravity. For 



example, the actions of the consistent supergravity truncations described in |32 
can be described by the action fl6"TT|). 

The equations of motion which follow from ( |6.1|) are 



1 



= ((dm-^A m )V^D^) 



19Vf ^ ■ dV ^\D^ + 



F z + 



dV 

dip* 



^d rn (sf=g~V F F mn ) = iqV^ (r(D n iP) - ip(D n ip)*) - C ^ e nm ™ s F mp F qs 



Rmn 2 ^9mn 



6g, 



mn ' mn 



T 

I in 



(6.3) 
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with 

Z o z z 



2 



+ ^ ((AnVOPnVO* + (D n ^)(D m ^y) . (6.4) 

Our strategy for computing the parity odd corrections to the entropy current and 
the parity odd transport coefficients closely follows that of |24], [26|, |2l], 116, |l(| . We 



first construct a stationary solution to the equations of motion ( |6.3|) in which neither 
the superfluid not the normal component are in motion. Being interested in the limit 
of small relative superfluid velocity, we then allow for the gauge field to have a small 
non vanishing spatial component. By boosting and rotating the latter solution, we 
find a generic bulk configuration involving six free parameters which are associated 
with the temperature, chemical potential, velocity of the normal component and the 
expectation value of the goldstone boson of the dual superfluid. We then promote these 
six parameters to be spacetime dependent fields. Once we do so the equations of motion 
(|6.3|) will no longer be satisfied. Therefore, we add corrections to the our previous 
solution so that the equations of motion (|6.3|) will be satisfied at least up to second 
order gradients of the six parameters. We can then compute the one point functions of 
the energy momentum tensor, charged current and gradient of the Goldstone boson in 
the boundary theory and from them read off the transport coefficients associated with 
the first order terms in a gradient expansion. 19 

6.1 A stationary solution 

The ansatz for a stationary solution to the equations of motion ( |6.3| ) — corresponding 
to a configuration in which the normal fluid and superfluid are both at rest — is 

ds 2 = -r 2 f(r)dt 2 + r 2 dx 2 + a(r)dtdr (6.5) 
^ = p (r)e^ (r) (6.6) 
A m = (A (r),0,0,0,A 5 (r)). (6.7) 

In this ansatz we have not completely fixed a gauge. In what follows we will present 
our analysis in the gauge invariant variables 

G m = A m - d m tp . (6.8) 

In this section we will not attempt to find an explicit solution for the unknown 
functions in (|6.5D, such a solution will not be needed in what follows below. We will 



19 The authors of Q recently argued that the transport coefficients associated with the parity odd 
sector can also be computed using a Kubo formula. 



- 51 - 



instead note certain abstract aspects of this solution that can be obtained without 
explicitly solving the equations. 



The equation for G5 is algebraic and is given by 

G. = -%. (6.9) 

The remaining equations of motion can be chosen to be a set of 4 linear differential 
equations in the variables /, a, p and Gq. The equations have the property that 
they involve only first derivatives (in the variable r) of / and a and second order 
derivatives (in the variable r) of p and Go- 
Assuming that the solution is a black hole, the fields a and p are non vanishing 
at Th, the black hole horizon, while / and Go have a simple zero at r^. The fact 
that Go vanishes at follows from regularity of G5 (defined in ( |B.9| )) near the 
horizon. The temperature, T, and entropy density s of the black hole are given 
by 

T = T } n " h \ s = ^4. (6.10) 

One linear combination of the 4 equations can be shown to express the constancy 
(in the radial direction) of the 'Noether' charge Q\ [BD], 



= (6.11) 
a a 

The equation in question asserts that 

Q[ = 

where a prime denotes a derivative in the r direction. Using ( |6.10 ) we find that 



Q 1 = sT. (6.12) 



• It is not difficult to work out a Graham Fefferman style solution of the equations 
of motion at large r. Putting in the physical constraint that the non normalizable 
mode dual to the scalar is turned off (i.e. that the operator dual to the scalar field 
is not sourced in the dual boundary theory) , we can use ( |6.11[) and the equations 
of motion for /, a, p and Go to show that near the asymptotically AdS boundary 
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(r — > oo) 

f = 1 + ^ + °(^) 

^ = l-^A|(0,)| 2 Ar- 2A + 0(r- 2A - 2 ) (6 13) 

p = r- A (C A |(0^)| + 0(r~ 2 )) 
Go = /i - n 2 q t r~ 2 + 0(r- 3 ) 

where Ca is a real number, O^, is the operator dual to if), and q t represents the 
total charge density as seen in the rest frame of the normal component; using the 
standard AdS/CFT prescription |36], |35| to convert bulk to boundary quantities 



we find that 

q t = -u^J^ 

with J M the boundary theory charged current. Using ( |6.11| ) together with 
we find 

sT + m = ¥ ^^^P (6-14) 
which is the Gibbs Duhem relation for a conformal superfluid (see (|3.6| )). 

6.2 Adding in a small uniform superfluid velocity 

As we have explained above, in this section we focus on the collinear limit, i.e. the limit 
in which the normal and superfluid velocities are equal. Following a long tradition we 
do, however, wish to allow the derivatives of the normal and superfluid velocities to be 
independent variables. Thus while we can set the normal and superfluid velocities equal 
at a given point, it would be inconsistent to do the same in the neighborhood of that 
point. In order to implement the fluid gravity map, even in our limited context, we need 
more general stationary background solutions than the ones described in the previous 
paragraph. In particular, since we are interested only in first derivative corrections, we 
need control over background solutions in which the fluid is at rest and the superfluid 
in motion to linear order in the superfluid velocity. In this section we describe the 
relevant solutions. 

A constant superfluid velocity is an 5*0(3) vector. The solution that describes 
infinitesimal superfluid motion (at first order) is a regular, normalizable, translationally 
invariant fluctuation about the static background described in the previous subsection. 
A vector fluctuation involves linear fluctuations of the vector modes gu and G,. We 
find it more convenient to work in terms of the variables g and 7 defined by 

Gi = -gdi<j> g ti = -r 2 7<%0 
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where di<p is the expectation value of the spatial component of the Goldstone boson in 
the dual field theory (i.e. d^4> is a function of the four boundary rather than the 5 bulk 
dimensions) . 

From the Einstein equations we obtain the following equation of motion for 7 at 
linear order 

' ir 5 + G' gr 3 V F Y 



a 



0. (6.15) 



Integrating fl6.15|) we obtain 



2 7V 5 + G' gr 3 V F 
2k z Q 2 = (6.16) 



with Q 2 an integration constant. The asymptotic behavior of 7 can be related to 
thermodynamic quantities in the solution via 

l = \{qt-q)r- i + 0{r^). (6.17) 

Integrating ( |6.16| ) near the boundary, it is not difficult to verify the form ( |6.17| ) and to 
check that 

Q 2 = q- (6.18) 

The equation of motion for g is obtained from the spatial components of the 
Maxwell equation and can also be written as a total derivative, 

/r3(gG °/ G ° )1/F + 2^,)' = 2k 2 /% . (6.19) 

Integrating it once we get 

fr 3 (gG' -g'G )V F , o , ^ _ ^ , - 



a 



lK l ^Q x = Q 3 + 2n 2 fQ 2 (6.20) 



with Q3 an integration constant. We will now argue that Q3 vanishes. The argument 
follows by demanding regularity at the event horizon. The fact that is an event 
horizon (and so that dr is a null one-form at r = r^) implies that 

7(n0 = 0- (6-21) 

Every term in (|6.20|) vanishes at r = r h except Q 3 , so we conclude that Q 3 = for 
T > 0. The near boundary expansion of g then takes the form 

g = 1 _ (qt-q)K 2 ^ 2 + 3 
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The second order linearized equations ( 5.15 ) and ( |6.19|) have four linearly indepen- 
dent solutions. One of these is 



g = 7 = constant (6.23) 

which corresponds to a deformation of the dual field theory metric. Another solution 
is 

9 = G 7 = tt/ (6.24) 

which corresponds to a boost of the stationary solution. Of the remaining two linearly 
independent solutions only one can be chosen to be regular at the horizon. In terms of 
the two first order integrated equations (|6.16| ) and ( |6.20| ), imposing Q 2 ^ precludes 



(|6.23|) and imposing Q 3 = removes the solution which diverges at the horizon. 



6.3 Boosting 



So far we have worked in a frame where uF" = (1, 0, 0, 0). As explained in |L6[ we can 
always boost to a frame where the metric and gauge field take the form 



(6.25) 



ds 2 = -r 2 f(r)u^u u dx^dx v + r 2 P^ u dx^dx v - lau^dx^dr + 2r 2 ^u^C, u )dx^dx h 
= Gu^dx^ + gC^dx^ + G 5 dr , 

with 

G = Gq 

lim A/3 = a p (6.26) 

ap the external gauge field and <fi the Goldstone boson. Our conventions here slightly 
differ from the main text where we've defined ( = —(, see equations flOp, ( |3.2j ). 

In our current holographic formulation, it is convenient to use a frame where the 
energy momentum tensor, charged current and Josephson condition are given by 



T^ = ( p + P + ftfuW + if v P + VK{,u u) + fUu + f ( 



/(;/ 

diss 



J" = + fC + JLs (6-27) 

(<9 M - a M ) = -11 + fidiss = -fi>T 

where 

<k = q + f», (6-28) 
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and our frame choice is u^T^ = and u^J^ iss = 0. 20 One can check that expanding 



(|6.27|) to linear order in and neglecting T^ ss , J% iss and jldiss, the boundary theory 
energy momentum tensor, charged current and chemical potential associated with the 
solution ( p.25| ) are obtained. 



6.4 The gradient expansion 

We now have in hand all stationary solutions that will be needed for our analysis. As 
explained at the beginning of this section, the next step in our analysis is to allow 
the thermodynamic variables to depend on the spacetime coordinates and look for 
corrections to the metric, 5g mn , gauge field, 6A m and scalar Sip, so that the equations 
of motion ( |6.3| ) are satisfied. Since we will be interested in the collinear limit, we will 
set ( = but keep its derivatives non zero. In this limit, the equations of motion for the 
corrections to the metric, gauge field and scalar naturally decompose themselves into 
scalar, vector and tensor equations under the 5*0(3) C SO (3, 1) symmetry retained by 
the background. We write the corrections to the metric in the form 



1 

3 ( 

- r 2 Sf u^dx»u v dx u + 2r 2 5V^ u v dx^dx u + r 2 5n^ u dx^dx v , (6.29) 

the corrections to the gauge invariant combination G m in the form 

5A m dx m = -5G Ufl dx^ + 5g^dx» + 5G 5 dr , (6.30) 

and corrections to the magnitude of the scalar as 5 p. 

The equation of motion for SG^ turns out to be an algebraic equation which can 
be solved. The remaining kinematic equations for SG, 5p, 5f and 5a are coupled. The 
vector sector contains two coupled equations involving 5V^ and Sg^ and the equation 
of motion for the tensor mode can be solved for since it is decoupled from the 



Sg mn dx m dx n = —2u^dx^r I u a d a u u + -d a u a u u ) dx v — 25a u^dx^dr 



20 Actually, in our computation we are using a frame where 

T>» = (p + P + ff)u»u v + rTP + 2/K(^) + — M* + TZs 

fix 

u M (<9 M - a M ) =-fi + ftdiss = -Mt ■ 

The dissipative corrections Jdiss, Tdi SS and fidiss in this frame are different from J diss, Tdi SS and Hdiss 
used throughout this work. Since in this section we will only be computing the transport coefficients 
in the vector sector, the difference between the two frames is immaterial. The reader is referred to 



16, [L0| for an extensive discussion. 
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rest of the equations. We have verified that the five constraint equations imply energy 
momentum conservation and current conservation in the dual field theory. We have 
also checked that the ratio of the shear viscosity to entropy density retains its universal 
value as has been alluded to in ]5l], Q [RJ. 



Since we are interested only in the parity odd sector, and since in the collinear limit 
this sector receives contributions only from the vector modes, we can focus entirely on 
the equations of motion for 5V^ and 5g^ which read 

pgg_gywy =sg , (ola) 

JMMlM.^)'^. (6 .31b) 

(Recall that G = G from ( |6.26| ).) Note that the homogeneous parts of these equations 
agree exactly with fl6.15| ) as expected. 

After integrating these once we can eliminate 5V in place of 5g and write the 
resulting differential equation in Sturm-Liouville form, 

K - — + tfraVjf? ) 5 gil = --£- + ^- / Si v (x)dx = S total . 



a \ a * vr j ^ G r 5 g J r ~» 

(6.32) 

In order to simplify the term linear in Sg^ we used the equation of motion for G. 
The limits of integration on the right hand side of fl6.32| ) have been determined by the 
boundary values of Sg^ and SV^ and the requirement that we are in the frame described 
in (pi. 

Since ( |6.32j ) is a linear second order differential equation we can use the method 
of Green's functions to solve it. Consider first the homogeneous version of Q6.32 ). We 
denote the two linearly independent solutions of the homogeneous equation by a and 
a such that 

a = 1 + 0(r- 2 ) a(r h ) = a h (6.33) 

a = r~ 2 + . . . d(rh) = oo . (6.34) 



We find that 



Qig + Q 2 G 



(6.35) 



Qi + Q2H 

The solution ( B.35 ) corresponds to a shift in the superfluid velocity described in the 
previous subsection. (One can use the method of variation of parameters to find a in 
integral form though we will not be needing it in our analysis.) The solution to 

d fr 3 fV F dg(r,?)\ frV 2 G' 2 o2 > 



dr V a dr J V a 



+ 2q 2 raV i ,p 2 Q{r, f) = 5(r - f) (6.36) 
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which vanishes at the asymptotically AdS boundary and also at the horizon is given by 

S{rtf = \ - (6-37) 

\—±a{r)a{r) r < r . 

Thus, if S' total (oo) = 0, a near boundary series expansion of 5g^ is given by 

Sg " = iy s r'"* + & js» (^|*r (*>) +■■■■ < 6 - 38 > 

Then, according to the standard AdS/CFT prescription, the dissipative corrections to 
the charged current in the boundary theory take the form 

J -" = ~h (f aST>dx ~ Ss. (^'"w) ) ■ (M9) 

Let us focus on parity odd contributions to S total . we find that 

S^ otal = Ack 2 (G - w)G'u, - 2cn 2 gG'B, + ( pa ^ r y m e s ven ) (6.40) 

where and were defined in Table [l| We reproduce their expression here for 
convenience 

^ _ X _^ P ° UvdpUa B » = -eT^F^ . (6.41) 

Hence 

= gi + C g^ (/"(^ + QiG)gG'drB» + jT 2(/x</ - G)(Q l9 + Q 2 G)G'druj t ^ 

, f parity even\ /g 42 ) 
V terms J 



Following the work of [26|] , the contribution of the parity odd terms to the entropy 
current come from the horizon value of 6V^, i.e., if we denote the entropy current by 



Jg then 



JS = *5V"(r h ) + ( pax 4 y m e s ven ) • (6.43) 

It remains to evaluate the parity odd contributions to SV 1 at the horizon. Integrating 
fl6.31b|) once we find that 

1 y -^-2n 2 Q 1 5V li = / S^ + SQ (6.44) 
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with 5Q an integration constant. Evaluating Q6.44 ) at the horizon and at the boundary, 
we find that 

jn [°°S 59 dr f \ 
XWVf*. \ — lL diss , Jrh I / parity even\ 

5V {rh) ~ OT + ~2Q^~ + V terms ) 



T 



JLs + f QH 2(H9 - G)GG'dr^ + J°° gGG'drB^ + (^^J en ) ■ 

(6.45) 



Restricting ourselves to parity odd contributions, the first term on the right hand side 
of the first line of ( 6.45|) represents the contribution of the canonical part of the entropy 
current </g canon - The second term on the right hand side of the first line of ( B.45| ) will 
give us the corrections to Jo 

O o canon 

Following the notation in ( 4.21| ), we denote 

R= jh (6 - 46) 

and write the parity odd transport coefficients and the corrections to the entropy current 
as in flOp 



_ z- , ,A» z- ra» _i_ ( parity even\ 
J diss - - ^bB + \v t J mg ) 

tA 4 I 1 tm \ st , ^ \ st , ,v \ ( parity even\ 

J s - -j; J diss + + o^uj + \ v te ^ mg J 
Then, (|6.45 ) and ( |6.42| ) imply 

POO 

k B = -c g 2 G' + R(G - gfi)gG'dr 



(6.47) 



(6.48a) 



POO 

k B = 2c\ (G - fig)gG' + R(G - fig) 2 G'dr (6.48b) 
2c ro ° 



POO 

/ (G — fig)GG'dr (6.48c) 

C f°° 

a B = -/ gGG'dr. (6.48d) 

1 Jrh 

While we can not solve the integrals in ( |6.48| ) explicitly, we find that the following 
relations are satisfied: 

1 CIJ? 

-<t u - fKT B = — (6.49a) 
1 _ cRu 3 

-k u - iik b + T (1 - fiR) o B = — . (6.49b) 
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(To obtain ( |6.49| ) from (|6.48|) we used G(oo) = fi and G{rh) = 0.) The first of these 
relations is precisely the universal prediction ( |5.27|) between entropy current coefficients 
predicted in the last section. Using ( |5.28|) and (|5.27|) we can also obtain holographic 
expressions for s$, a\, as and a\o 



s 9 = 



g(AG - ng)G'dr 



g'G'dr 



aio = J g(G-fig)G'dr 



(6.50a) 
(6.50b) 

(6.50c) 
(6.50d) 



To study the values of k w , Kb, and <jb at the phase transition we set g = 1. 
Then, all the integrals in ( 6.48|) can be carried out explicitly. We find that the resulting 
expressions exactly match the zero superfluid density values of k u , kg, cr w and as 
studied in |2T], ^ ^DJ , implying that these coefficients are continuous across the phase 
transition. 



7. Discussion 

In this paper we have described a framework for describing superfluid hydrodynamics 
at first order in the derivative expansion. We have determined the most general form of 
the hydrodynamical equations that are consistent with Lorentz invariance, the Onsager 
principle, and the requirement that the second law of thermodynamics apply in every 
conceivable situation. We begin this section by summarizing our results, and then turn 
to a discussion of possible applications and extensions. 

We have found that the most general equations of Lorentz invariant but non-parity 
conserving superfluids requires the specification of twenty parameters. Fourteen of 
these are associated with the parity even sector and were described already in JTD[ , 



generalizing an earlier 13 parameter framework spelt out in JTB|, |T^|). In this work 
we have found that six more parameters need to be specified in the parity odd sector 
of the theory. Five of these six parameters are non dissipative; they drop off from 
the expression of the divergence of the entropy current, and are unconstrained by 
inequalities. Three of these five parameters (fj, R n and £22 i n the language of the main 
text) are relatively simple. They are the transport coefficients which are associated 
with the pseudo-tensor mode 
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and the pseudo vector modes 



Vi = e^uX P [d£-^\ (7.1b) 



T T 

V 2 = e^u v ( p ( a <j aa . (7.1c 



Note that fj closely resembles the hall viscosity in 2 + 1 dimensional theories There, 
when parity is violated it is possible to have a transport coefficient rj H associated with 
the tensor mode 

e^ p u u a p a + e aup u u a p p . (7.2) 

The hall viscosity is known to be associated with the ground state's intrinsic angular 
momentum. It would be interesting to see whether fj is also associated with similar 
quantities. 

The remaining two nondissipative constitutive parameters, erg and a\Q, multiply 
relatively complicated expressions in the constitutive relations. These two parameters 
determine the eight transport coefficients associated with the pseudo scalars u;-£ and 
B-£ which appear in the constitutive relations, and the transport coefficients of the 
pseudo vectors *d tl T, *d fl ^ J ^d^^, w' 1 and B p . It is possible that these expressions 
admit significant simplifications when expressed in terms of variables that might turn 
out to be more natural than those adopted in this paper, as and <Tio also parameterize 
deviations of the entropy current away from its canonical form. 

The sixth parity odd parameter is dissipative since it enters into the formula for 
entropy production. In the main text we have denoted this parameter by k 2 i (and also 
k\2 = — «2i). It is associated with the pseudo vector modes (|7.1b|) and ( |7.1c|) . While 



undetermined, it, and the transport coefficients Kij i,j = 1,2 associated with the parity 
even sector, need to satisfy the inequality 



12 



Note that this equation mixes parity odd and parity even coefficients. 

The triangle anomaly, if it exists, enters the expression for the pseudo scalar and 
pseudo vector coefficients whose value depends on a 8 and a w . It is interesting to 
compare this result with that of [^D] which studied (non-superfluid) hydrodynamics in 
the presence of triangle anomalies. There, the only addition to the parity-odd sector are 
the two pseudo vectors u p and B p whose transport coefficients are completely fixed in 
terms of the anomaly and a single integration constants. When considering superfluids 
the functional dependence of the transport coefficients of u p and B p on T, /x and ( is 
completely arbitrary. This novel feature opens the possibility for interesting physical 
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phenomenon in a parity violating theory, even in the absence of an anomaly which is 
an intrinsically relativistic effect. 

The most immediate application of our work would be to the modelling of the long 
distance behaviors of real world superfluids (or superconductors) whose hydrodynamics 
violates parity. We do not yet know of any candidate experimental systems of this 
sort. We do note, however, that non centro symmetric superconductors live on parity 
violating lattices. It would be exciting to investigate whether parity violation in these 
(or analogous systems) could lead to experimentally observable nonzero values for any 
of the 6 parameters described above. 

In this paper we have worked out the general theory of superfluidity in 4 dimensions. 
Several aspects of our analysis dependend on the existence of an e tensor with 4 indices. 
It seems likely that the theory of superfluidity in 5 and higher dimensions will differ 
in qualitative aspects from the theory we have worked out. It would be interesting to 
flesh this out. 

We have already verified certain aspects of our general construction of the equa- 
tions of superfluid hydrodynamics using explicit computations within the AdS/CFT 
framework. However all the computations reported in this paper work in the collinear 
limit of vanishing superfluid velocity. It would be instructive to demonstrate that all 
the numerous other relations, implied by our work, between transport coefficients that 
do not preserve SO (3) invariance, are also borne out by AdS/CFT computations. In 
this context we point out that the generality of the expected results suggests that their 
derivation may also be carried out for a generic holographic superfluid (perhaps along 
the lines of Section ^) and should not involve the details of the corresponding dynamical 
systems. 

An interesting feature of our AdS/CFT analysis is that in the collinear limit all 
the transport coefficients in the parity odd sector of the theory are continuous at the 
phase transition (i.e., they match their value in the uncondensed phase where / = 0). 
This is similar to the behavior of the shear viscosity r\ and diffusion coefficient k of the 
parity even sector, but differs from the divergent behavior of the transport coefficient 



associated with the scalar d^f^) [16|. It would be interesting to understand this 



result in terms of the theory of dynamical critical exponents [[54 

It would also be interesting to work out the properties of stationary superfluid 
flows (like rotating superfluids on an S 3 ) and to investigate how the properties of such 
configurations are affected by the parity odd non dissipative terms that can be 'turned 
on' in such configurations. Such configurations would be dual, under the AdS/CFT 
correspondence, to rotating hairy black holes in global AdS space. 

To end this paper let us highlight a structural aspect of our analysis that we find 
quite remarkable. For this purpose it is sufficient to focus on the case of parity invariant 
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superfluids. The total number of constitutive coefficients allowed, merely on symmetry 
grounds, is 50. In this paper we have shown that the requirement that the equations 
that follow from these constitutive relations are consistent with positivity of divergence 
of any entropy current (plus the Onsager relations) cuts down the number of constitutive 
coefficients from 50 to 14. In other words, the second law of thermodynamics gives very 
powerful and precise constraints on dynamical equations. Given the duality between 
fluid dynamics and gravity, it is natural to wonder whether similar results may also 
be true for extensions of the theory of gravity. We leave detailed investigation of this 
exciting possibility to future work. 
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A. The linear independence of data for the parity even super- 



fluid 



A.l The linear independence of first order terms 



A. 1.1 The vector sector 



The equations of motion in the vector sector are 



(A.l) 
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A basis of ten one derivative vectors (before using the equations of motion) was listed 
in table [| It is given by 



(A.2) 



P^(u-d)up, P^(u-d)^, M{^8)u Pi P^dp 
ptfdpT, P^dp (J£\ , P^Cd0U Vi P^Ep, P^F^C- 



The quantities in ( |A.2j) are not all on-shell inequivalent as they are constrained by the 
relations ( A.l ) . In this subsection we will argue that it is consistent to choose the seven 
vectors listed in the third column of Table ^ as independent vector data. That is, we will 
show that it is possible to use the equations ( |A. 1| ) to solve for P^dpT, P^dp (|;) , P^^dpUi 
in terms of 

P^(u-d)up, P*{wd)fa P*(Z-d)up, P*{t-d)fa P^df, (£) , 

If we rewrite the equations of motion in ( |A.1| ) in terms of the quantities in (|A.2 ) 
we find 

P^ ((P + p){u-d)up + d T PdpT + df>d p £ + + =P* (qE, - fF^C) 

(( u -d)& - Td p £ - ^dpT + fdpuu) =P^E fi 

f + + ^dpTj =P« , FfW. 

(A.4) 

It is possible to use ( A.4 ) to solve for the scalars listed in ( |A.3| ) if and only if the 3x3 
matrix of the three vectors P^dpT, P^dp (|;) , P^^ v dpu v in the three equations 
( |A.4| ) has nonzero determinant. This 3x3 matrix is given by 

d T P % T P ' 
M (v) = | —fi/T 1 | (A.5) 



e/T Ti 0, 



and its determinant is given by 



Det (M {v) ) = f f^P - Td T P^j . (A.6) 

It is nonzero for a generic functional form for P(T, //,£). We conclude that the vectors 
( A.3|) form a basis for onshell independent one derivative vectors. 
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A. 1.2 The scalar sector 

The equations of motion in the scalar sector are given by 

i v d^ = qE-t 

(A.7) 

u v ^{d^ v -d v Q = E-i. 



u v d,T^ = fE-i 
= cE-B 



A basis of 11 one derivative vectors (before using the equations of motion) was listed in 
Table ^. We denote them by {C^\ S\ } for the first set of on shell independent scalars 
and {£-f \ 5$ } for the second set. Here j runs from 1 to 4 and % runs from 1 to 7. We 
have used the notation in Table £|. The new quantities Cj are defined as follows 

4 6) = £-0£i, (A.8) 
4 6) = £^-<9<v 



£ (a) 


= U 




r (a) 
'-2 


= U 




r (a) 
^3 


= u 




r (a) 
'-A 


= ^u-du 



The quantities defined in ( |A.8|) are the dependent data for the two choices of bases 
among the on-shell inequivalent quantities. These quantities are to be determined by 
the equation of motion ( |A.7| ) in terms of the dependent quantities 5*,. Note that the 
sets {C[ a \ } and {Cf \ S^} are different partitioning of the same set of quantities. 

The equation of motion in ( |A.7[ ) expressed in terms of the quantities in ( |A.8| ) has 
the form 

E(eS a) ) P sf } + E(4 a) ) P q = °- 

T 7 (A.9) 

i=l i=l 

In the equations above the index p runs from 1 to 4 denoting the 4 equations in ( |A.7| ) . 
Again note that both the equations in ( |A.9| ) refer to the same set of equations. We find 
it convenient to define the new set of quantities 

A = ~t 2 {u 2 -x 2 y B = ~T 2 {^-ey c = ~t(z/ 2 -e)' ^ A ' 10 ^ 

so that the projector 

= jf + Au"u v + BCC + C (£V + vP?) • (A. 11) 
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The coefficients in (|A.9|) are given by 



and 



(a) 


)l 


= "(P + p), 


(a) 
e 3 


)l 


= B(P + p)+Cfp-f, 


(a) 


)l 


= pd u f-CTfp + fT, 


(a) 
e 7 


)l 


= -f 


(a) 
e i 


)2 


= (p + p)h, 


(a) 

er 


)l 


= -(pB(p + P) + CTf£ 


_(a) 
e 5 


\ 

J2 




(a) 




— — n 
Hi 


(a) 


)s 


= q, 




)3 


= -(Bq + Cf), 


(a) 


)s 


= dJ + CTf, 


(a) 
e 7 


)3 


= 0, 


(a) 


)4 


= 0, 


(a) 
e 3 


)4 


= -1, 


(a) 


)4 


= T, 


(a) 
e 7 


)4 


= -1. 



(e? 
(e? 



fa 



fa 



fa 



fa 



B^Tfp + pd x f, 
pd T f + 2fu, 



(A.12) 



«9 r P - ed T f - 2f X 2 T 



d x P-C 2 d x f-UT-BeTf, 



I 
T 



-d x f - BUT 
-d T f - t 



) 4 = o, 
) 4 = o 

)e = v, 



(4 a) )i 


= Cfp£T-d x p, 


(4 a) 


)i 


= -(AfpT + d uP ), 


(4 a) )i 


= ~d T p, 


(4 a) 


)i 


= Afp + C(P + p), 


(4 a) ) 2 


= pd x (p+ P )-cfeT, 


(4 a) 


) 2 


= pd u (P + p)+TAef, 


(4 a) ) 2 


= pd T (p + p), 


(4 a) 


h 


= (P + p)-Cp(P + p) 


(4 a) ) 3 


= 8 x q - CUT, 


(4 a) 


h 


= d u q + ATf, 


(4 a) ) 3 


= d T q, 


(4 a) 


h 


= -(Cq + Af) 


(4 a) ) 4 


= U, 


(4 a) 


h 


= 0, 


(4 a) ) 4 


= X 2 T, 


(4 a) 


)4 


= 0. 



(A.13) 



Afe 
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The other set of coefficients, with index (b), can be read from ( A.12|) and ( |A.13| ) using 



(6), 



(b)\ 
2 Ji 



<4 



(o) V- ( e {b) )- 



(4 



(*>)\ 
3 /« 



"6 A 



"3 A 



(4 a) ) 



4 A) 



(A.14) 



We can express all the derivatives in ( |A.12 ) and ( A.13|) as derivatives of a single func- 
tion, say, the pressure. Thermodynamic relations that enable us to do so are 



q = j, 9 h/tP ; / = J^d i/T P 
We make the following observations: 



-P + Td T P-j;d^P. (A.15) 



a) We can use the equations of motion ( |A.9j ) to solve for the 4 scalars ^u-Vu^ u-dYn 
(i = 1, . . . , 3) in terms of the 7 independent scalars in the 3rd column of the first 
row of Table |[ This is possible if and only if the 4x4 matrix of coefficients of 
the four quantities in the first equation in (|A.9|) has nonzero determinant. This 
matrix is given by 



(a) 



ij 



(A.16) 



b) We can use the equations of motion ( |A.9|) to solve for the quantities ^^-du^, 
£-<9Ej in terms of the 7 independent scalars in the 3rd column of the second row 
of Table |3|. This is possible if and only if the 4x4 matrix of coefficients of the 4 
quantities in the second equation in ( |A.9| ) has nonzero determinant. This matrix 
is given 



M. 



u 



(4 



IJ ■ 



(A.17) 



The relations ( |A.15| ) allow us to express the matrices ( |A.16| ) and ( |A.17j ) in terms of 
the pressure. Using several reasonable equations of state we have used Mathematica to 
verify that the determinant of the matrices in ( [A. 16|) and ( |A.17|) is generically non-zero. 



A. 2 The linear independence of the second order terms 

A list of second order scalar data, the second order equations of motion and a choice of 
second order independent scalar data can be found in Table The second order scalar 
equations that follows from the first order vector equations are 



V a ( P a ^ (d^ u - d^u)) = V Q (P^F^ 



(A.18) 
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The two derivative terms in these equations take the form 



(A.19) 
(A.20) 
(A.21) 



The quantities P^V^d^T, P^V ^dv^, P^^V ^V^Ufs can be solved for using equations 
([A.19 ), ( A.20Q , and ( A.21 ). Note that these two derivative scalar quantities do not 
appear in any other equations of motion. We can then use the remaining 8 equations 
of motion to solve for the other 8 dependent data, 

u'VV^Ei, uW^V^up, eCZ^tNvUp, eCV^i (A.22) 

where % runs from 1 to 3. The reaming 8 two derivative scalar equation of motion are 







= uN p 


(- V) , 


(A.23a) 










(A.23b) 




u p V fi (V M J") 


= vPVp 




(A.23c) 










(A.23d) 






= eve 




(A.23e) 






= e /3 v /3 




(A.23f) 










(A.23g) 










(A.23h) 



The matrix of coefficients of the terms in QA.22 ) as they appear in the equation of 
motion ( |A.23| ) may be expressed as 



Nii 



M$ 




(A.24) 



where the rows represent the ordered equations in ( |A.23| ) and the columns represent 
the ordered quantities in ( A.22Q . It follows that 



Det [Nij] = Det [M^f } Det [M % 



r(&)i 

-ij J 



(A.25) 



In the previous section we concluded that both Det[M^] and Det[M i i0 ' 1 ] are generically 



XI 



non-zero. Therefore we can infer that Det[iVy] is also generically non-zero. 
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In order to understand the structure of the matrix iV we note that the first four 
equations in (|A.23|) are generated by the action of u- d on the first equation in ( |A.9|) . 
We then find that u-d acting on generates all the independent second order data as 
presented in Table |5]. Likewise, the action of u- d on the C\ generates the four terms 
u^u v ^V v up, u^u y V \td v Y,i {i = 1,...,3). In fact these dependent two derivative 
terms appear only in equations ( |A.23a| ), ( |A.23cj ), ( |A.23d| ), ( |A.23e[ ) and is not there in 
the rest of the four equations in (|A.23|) . 

Similarly, we can think of the equations ( |A.23f| ), ( |A.23g| ), ( |A.23h| ), ( |A.23h| ) as being 
obtained by the action of £• d on the second equation in ( A.9 ). Also here the terms 
V M V^, f^" V M d„£< (which constitutes the 4 remaining second order quantities 
which are determined by the equation of motion) are generated by £ • d acting on the 
Cf' terms. These dependent four second order quantities do not appear in the first four 
equations in (|A.23 ). This structure justifies the block diagonal form of the coefficient 
matrix in ( |A.24|) . 



B. Derivation of thermodynamic identities for normal fluids 



In this appendix we present a derivation of the equations of motion ( |2.22| ). Our starting 
point is the observation that the scalar components of the equations of conservation of 
the stress tensor and current reduce, at first order to 



u v V„T^ = u u F^J u = 
(ud)p + (p + P)Q = 



(B.l) 



= (u-d)q + qQ = 



(B.2) 



where 9 = (V-u) 



Derivation of the first equation in ( p.22| ) 
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Subtracting [ eq. ( p.2| )1 from [|^x eq. (|B. 1|) ] we find the following relations: 



- {u .d)p--{u-d) q 



dq dp 
dv q dv 



dq dp 
pv^aT 


dp dq 
~ dv~dT_ 


(u-d)T = - 


(P+P) 


dq dP 


dP dq 


(u-d)T=- 


'dP~ 




~dv~~&f 


~ ~dUdr 


dp 


g 



e 



, dq dp 



e 



p + P\ dq 



T 

u-d)T 
T 



dv q dT 



(u-d)T=- 



~dP~ 




dp 


<? - 



{p+p) dv-- q dv- 



e 



dP 
dp 



e 



In the third line we have used two identities: 



dP 




dp 




dP 




dP 




dq 


dp 


<? 


dT 


u 


dT 


V 


dq 


p 


dT 


~dP~ 




dp 




~dP~ 




~dP~ 




dq 


dp 


q 


dv 


V 


dv 


V 


. &q . 


p 


dv 



(B-3) 



(B.4) 



In the last line of ( |B.3|) we have used the following three thermodynamic identities. 



dP 
~&T 
dP 
dv 
dp 
dv 



p + P 
T 

Tq 



(B.5) 



T 



d 2 P _dP__ 2 d_ 
dTdv ~ ~dv~ ~ dT 



1 dP 
T~dv~ 



T 



dq 
df 



Derivation of the second equation in (|2.22|) 

Adding[ gx eq.flB^)] and [gx eq.([BlD1 we find 



dP dP 

— (u-d)p + —(u-d) q 



n ,dP dP 



e 



(u-d)P = - 



nN dP dP 







(p + P)^^ r - + Tq(u-d)v 



{u-d)v 



T 
1 

f 



nS dP dP~ 
ip + P) ^ + q ^. 



(B.6) 







~dP~ 
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In the third line of (IB. 61) we have used the first law 



dP=(S + uq) dT + Tqdv = (p + P) — + Tqdv 
p + P = T(S + uq) 



(B.7) 



and in the last line of ( |B.6| ) we have used the first equation in ( 2.22| ). 



Derivation of the third equation in Q2.22] ) 

The third equation in (2.22) follows from the vector component of the stress tensor 
conservation equation. 



PtfVJT* = P, e F ea J a = qE, 
^P°d e P +(p + P)(u.V)tv = qE^ 

Q 



Vi, 



(B.8) 



In the last line of ( |B.8| ) we have used the first law as written in ( |B.7| ). 



C. Details relating to parity violating superfluids 

In this Appendix we provide several computational details relating to section [§ 

C.l All equations of motion for ideal superfluid 

The equations of motion for a superfluid are listed in table [| There are four scalar 
equations, one pseudo scalar equation and 3 vector equations. In this subsection, using 
thermodynamics, we simplify these equations so that they can be easily used to solve 
for the dependent fluid data in terms of the independent data. 
First we note the identity 

^-£v^-(£-v)£ M =i^r- 



In the second line we have contracted both sides of the equation with £„. Using flC.ip 
and the first law we can simplify the stress tensor conservation equation projected in 



the direction perpendicular to m m : 
P»V u T ud - P^F du J u = 



=^V„ 



(p + p)u»u» + Prr + /re 



>p^ e [v,p + /(£• V)6 + /F^ei + (p + VK + C M v e = 

>P^ [V,P - fZVet] + (p + P)(ii-V)tt" + C M V, = gP^ 



(C.2) 



(p + P)^ + TgV e z/ 



+ (p + P)(«.V)«" + CV^ (/£") = gP A 



V.P 
T 



(w-V)?/ 



p + P 



p + p 



In the second step of (|C.2|) we have used the identity ( |(J. 1| ) and in the third step of 
(C.2) we used the first law in the form: 

dP = (p + P)— + Tq du + fi d£. 

Next, using ( |C.2j ) we can simplify the equation relating the curl of the phase velocity 
to the field strength (ie. C^ u = — — F^ u = 0). Consider the case where 
one of the two free indices of this equation is projected in the direction of w M . 



u p (V^-V^-P, 







=W [V M (-p u p + Cp) - V^-p + Q] = P M 

[V M p + u»{u- V)p] + p(n- V)^ - («• V)C„ - C e V M M e = P M 
=> [PfteH + p(«- V)«J - (u- V)C M - C'V M n e = P^ 
=^V fl i/ + i/ [P M e V,T + r(wV)wJ - (w V)C„ - C V^u e = E„ 



(C.3) 



>( M .V)C M = -T 1 



pg 



p + P 



In the last step we have used 

Next consider the case where both the indices of are projected in the direction 
perpendicular to u^. This can be simply analyzed by contracting the two free indices 
of C mj with e^ A<T w„ 



e^u v (V A ^ - V^a - Pa.) = 
^2e^u u V x ^ - e^u u F Xa = 

- 2pe^u u V x u a + 2e^u u V x ( a - e^u u F Xa = 

^ ^ = — + pu/\ 

/Li 



(C.4) 
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where f2 M , uj^ and are defined as 
1 



2 2 

If we project flC.2| ), (|C.3|) and (|C.4|) in the direction perpendicular to e , it will be a 
rewriting of the three vector equations as listed in Table |3|. Contracting the free index 
with e we get two of the scalar equations and one pseudo scalar equation from Table 
H rewritten in our basis. 21 

Next we rewrite the two remaining scalar equations of table |3|. First, using these 
two remaining scalar equations, we will show that the entropy current is conserved in 
equilibrium. It is most useful to use the basis in Table 







( P + p)uv + Prf v + fee 



+ fu v F v ^„ = 



- (u-V)p - (p + P)(V-u) + fiVeifC) + fu v [(£• V)e + F^Q = 

- [(u-V)p + /£(«• V)£] - (P + P)(V-u) + fxVe(fe) = 

- [T(u-V)s + fi(u-V)q] - (Ts + fiq)(V-u) + pN 9 {f£ 6 ) = 

- TV, (su^) - ^V M (gu" - fe) = 

In the fourth line we have used the identity ( |C.1| ). In fifth line we have used 

dp = T ds + p dq — f^ d^ 

and 

p + P = Ts + pg 
In the last line we have used the fact that 

V M J^ = V^qu" - fe) = C(Two derivatives) ~ 

Using V m (sm m ) = and V M (gw M — fe) = we find 

v ll [fe} = s(u-v) [| 



(C.5) 



(C.6) 



21 The relevant equations are 
1. CmV^T"" - C^""^ = 



2. u^CC^ = 

3. £^ Act C^Ca ct = 
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To summarize, we list all the equations in simplified form. The four scalar and one 
pseudo scalar equation are given by 



(V-u) = = 



(u-d)s 



v e (/e e ) 

(P + P) 



K = 



[B 1 (u-d) X + B 2 (u-d)u + B 3 (u-d)T] 
s(u-d) (§) 



[K 1 (u-d) X + K 2 (u-d)u + K 3 (u-d)T] 



K ' X ir + ( (> (u-V)u e = RT(V 1 -()-( 2 K 



T 



T 2 + A 2 (u-d)v + A 3 (u-d)T] - T(l - pP)(y r C) 



The three vector equations are given by 



T 



+ (u-vk 



F 7 "^ = P CTM 



RTP^V* 



+ pw CT 



P^ (u . V )c = p£ [-T(i - pP)\^ - cv^J 



where we have defined 



0l 

r 



X 



' 



K 



p+p ' 



p 



p + p 



Pi = 
= 

Ai = 

K = 



~pog(*)], P 2 = -|;[log( S )], Bs = -^pog(«)] 



9 



P + Pdx 
-1-^(1-^) 



5 



p + Pdv 



d 



p + PdT 



d (q 



- pP) 



- pP) 



#r Is J 



(!) 



+ 



X_ds_ 
3s dx 



+ 



X_ds_ 
3s dv 



X f ds 3s 
3s V <9T T 



^ - P a V r 



P 

LTJ 



(C.7) 



(C.8) 



(C9) 
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C.2 Showing the linear independence, of the first derivative scalar data 

Using the scalar equations we can solve for four dependent scalars and one dependent 
pseudo scalar in terms of the independent ones as they appear in the list of Sf in table 
^| We choose the five dependent scalars to be: 



e, k, c-<K, O-vk, n-(. 



(CIO) 



Note that the last term is the pseudo scalar. The dependence of the scalars in (|C.10|) 
on the iS 4 c 's from table |6] takes the following form 



= B 3 S 2 + B 2 S 3 + B\S1 

K = K 3 S C 2 + K 2 S C 3 + K^l 

(■a-( = A 3 S C 2 + A 2 S C 3 + A 1 S C A - T(l - iiR)S{ 

C{u-V)u v = RTSl + C [K 3 S C 2 + K 2 S C 3 + K^t] - %S\ 



(C.ll) 



(B-0 



where Ki and are defined in equation (|C.9| ). Using the first three of the above 
equations one can form the A s matrix in (|4.20|). 



C.3 Showing the linear independence, of the first derivative vector data 

We wish to argue that the vectors V c form a set of independent vectors. To do so we 
solve for the three vectors that do not appear in the list of in terms of using 
the three vector equations of motion. The three vectors that are not part of the set 
{Vf } are 

P^{u-V)u e P^iu-VXe P 6 >*{(-V)(e 
From the first vector equation ( |C.2|) we can solve for P etA -{u-V)ue ) 



(C.12) 



P^(u-V)u e = P 0fl + RTVie ) 



(C.13) 



1 T 



This equation has been used to form the matrix A v in (|4.20| ). 

Taking the star of the second vector equation (|C.4| ) we can solve for P dfJ, ((.'V)Ce 



p^(C-v)o = P etl 



T 2 

K a {d a u e - d e u a ) + — dex + T X d e T - F ex ( x 



2 



(C.14) 
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From the last vector equation (|C.3|) we can solve for P dfl (u-'V)(e 



P^iu-VXe = P e » -T(l - fj,R)V lg - CVofl - ^-{d e u a - d a u e ) 
= -T(l - (iR)V? - - C 2 V 6 CM . 



(C.15) 



C.4 Showing the linear independence of the two derivative parity odd scalar 
data 

Here we solve for the five two derivative parity odd scalars using the five parity odd 
two derivative equations of motion as listed in Table [7]. The solution will be presented 
in terms of the three independent two derivative pseudo scalars as listed in the last 
column of Table [7| and squares of single derivative terms. The two derivative pseudo 
scalars that we shall solve for, are the following 



e Xa (u-V)V x u a 
e A >-V)V A C 



e A,7 (C-V)V A C (J 



In the solution below we shall write the dependence on the two derivative terms explic- 
itly but will not write the terms which are squares of a single derivative. 



Ex = [e^ A XCJ VA(V e T a e - F a9 f) = 



Products of one derivatives 



( W .VK+y-^-T^) 



Products of one derivatives 



[e^ Xa u^ u ] {u e V 9 V x u a - RV\E a ) = Products of one derivatives 



(C.16) 



Using E% we can solve e X(7 {u-V)V \u a in terms of e A<7 V \E a . 



E 2 = [e^Q V u F Xa = 

[e^CtUu] (u-V)F Xa + 2 [e^CM u 9 V u F er7 = 
=>tv(u- V) [e^ Xa u u F Xa ] + 2 [e^C^ux] V u (u e F 9a ) = Products of one derivatives 
=>> [e^^u^ux] V x E a - C,^u v V u B^ = Products of one derivatives 

(C.17) 
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Using E2 we can solve ( fJ, u u 'V u B fl in terms of e Vj£ ff . 

£ 3 = [e^ A< %] V u F Xa = 

=^V M \e^ vXa u v Fxa\ = Products of one derivatives 
^V^-B^ 1 = Products of one derivatives 

^P^VpBv + ^-V^B V = Product of one derivatives 

Using E3 we can solve P^Vfj,B v in terms of ^V^. 
E , = [e^ A XC] mVaCW = 

^V A ( [e^XCj [{u-V)Ca + T{1 - fiR)V a + ( e V a u e ] J 



Product of one derivatives 



c 6 



{u- V)V A Ca + (1 - nR)V x E a + ^-[V A , V> 



= Product of one derivatives 
Using E4 we can solve e Acr (-u- V) VaCt in terms of e Xa V\E a . 

c \iv\a 



E 5 = [e^ A XC] (C-V)C Act = 



(C-V) VaC. - M(C-V) Va« ct 
= Product of one derivatives 



(C.18) 



(C.19) 



(C.20) 



Using i? 5 we can solve e Ao "(C.V)VAC<r in terms of e Xa ((.V)V \u a and ^^"W ^B u . Thus, 
we can choose C^Va-E^, e Ao "(C-V) Vam ct and ( fJ, ( u V^B u as the three independent two 
derivative parity odd scalar data. 22 

C.5 Relating the two different entropy current 

The relation between the variables in ( f4.3|) and ( (4.9|) can be obtained by using 
1 

— < 

2 



; e» uXa ^F Xa = -fiB" + (5<K - TV? ~ TV? 



l -^ x °V u [Tu x Q=T 



B^ 

T(^ + — 



(C.21) 



22 Here 
and 
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Explicitly, we find 



Si 


= Ta 1 




(C.22) 




= Ta 8 




(C.23) 


~s\ 


1 H l 
= a 3 - —(7i + — a 9 

s s 




(C.24) 


~b 

s 2 


T ii 


T 


(C.25) 


Vl 


RT 2 „ 
z 


(C.26) 




= « 2 




(C.27) 


v 3 


= Ta 3 + Td T (Ti 




(C.28) 


V 4 


= Ta 4 + Td v a x - Ta 8 




(C.29) 


V5 


= Ta 5 + Td x a x 




(C.30) 


v 6 


= —TfMTi + T 2 (Jq 




(C.31) 


v 7 


T 




(C.32) 



C.6 Computation of the divergence of entropy current 

V u Jq can be calculated term by term. 

o new *J 

V M (e^V u [Ta lUx C]) = (C.33) 



V M [TV 9 ^] = -2T"a 9 (u-{)K + 2KT 3 a 9 (u-V) 



+ T 2 



— ( W -VT) + -^-(w-Vi/) + ^Wx) 



(C.34) 



V^Tcno^] = -Ta w (B.()K + RT 3 a 10 (B-V) - 2T 2 a 10 (u-V) 



2TV 10 (w-Vz/) +T 



^( S .VT) + ^(RV,) + ^(5.Vx) 



(C.35) 



V M ( ( 7 3 V; 



'"^ = -2<7 3 (o;.C)( W .V)r- ^(Vf-VT) - ^(y^.VT) 

ov ox 



+ a 3 [i?T(^-V 3 CM ) + 2/i(wVT) + (S.VT)] (C.36) 
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V„(Ta 4 V 1 



Cf l \ 



-2Ta 4 (w-C)(M-V)z/ + T 



r; (v;".vv) • ^(V;".v\) 



+ Ta 4 [i?T(\/-V 4 CM ) + 2/i(wVi/) + (B.Vu)] (C.37) 



V M (To- 5 Vf ) = -27V 5 (u;-C)(u-V)x + T 



0(Tb {V?.VT)-^{V?.V X ) 



+ Ta 5 [RT{V-VT) + MwVx) + (B.V X )] (C38) 



V M (ye^e^A.) = (B-CXu- V)<7 8 - M^- V)^s - 27V 8 [(5- V) + 5- W)] 



+ T 



<9T 



9(7: 



da* 



<9z/ 



T 



<9x 

^(vr-vT) + ^(vr-v X ) 



In deriving these expressions we have used the following identities. 



RT 1 



-7 + — 



C.7 Frame transformation formula 

Under the frame transformation 



T' — T — ST 



v 



v — 5v 

the stress tensor, current and /idiss transform in the following way. 



(C.39) 



(C.40) 



(C.41) 



TTss ~ TZs = STZs = 6{p + P) + 5P rT + 6f + (p + P)(u»5u» + 
Jts ~ JLs = ^ J Ls = 5qu»- 5f e + 



Pdis 



Hdiss = Sf^diss = 5(j,-£-5u 



(C.42) 
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where 5 A for any scalar function A denotes 

r a OA dA 

= df 6T + to 5 " 

Using these expressions one can deduce how Sj transform under a change of frame 

5s x = 25P 

6s 2 = 5P + ( 2 6f 

5s 3 = 5p + p, 2 5f 

6s A = ^ 2 6f-(p + P)(6u-0 (C.43) 

5s 5 = -5q - fi5f 

5s 6 = -( 2 5f + q(5u-0 

5s 7 = —5fi + (5u •£) 

and also how the Vj transform under a change of frame. 

5^ = -(p + P)(P^5u u ) 

= (C.44) 
5v£ = q(P^5u u ) 

where '5' denotes a change under field redefinition. 

Using these formula one can easily form the frame invariant combinations in the 
vector sector. It turns out that is frame invariant by itself .The other frame invariant 
combination is proportional to 



q 



p + p 



+ v£ = R vf 



These are exactly the combinations that appear in the left hand side of equation ( |4.26 ). 



Checking the frame invariance of the combinations appearing in the scalar sector is 
more involved and require many thermodynamic identities. Here we have used Mathe- 
matica (version 7) to impose all these identities by expressing p, q and / in terms of a 
single function P(T, u, £) and its derivatives. Then we explicitly checked that the four 
combinations appearing in the left hand side of equation ( |4.35| ) are invariant under the 
transformations given by ( C.43|) . 
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